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EVALUATION 


This  contractual  effort  resulted  in  the  development  of  a geometrical 
theory  of  diffraction  (GTD)  computer  code  for  circular  parabolic  reflectors 
with  feed  and  strut  scattering.  This  computer  code  analytically  predicts 
the  far  field  pattern  of  the  parabolic  reflector  including  the  scattering 
from  reflector  edges,  feed  supports,  and  the  feed  structure. 

The  results  of  this  effort  provide  a means  of  simulating  a parabolic 
reflector  with  obstructing  feeds  and  supporting  structures  and  predicting 
the  far  outside  lobe  responses.  This  computer  code  can  more  accurately 
simulate  the  parabolic  reflector  environment  than  previously  accomplished 
in  the  past. 

This  fits  into  the  RADC  Technology  Plan  (TPO  I-B)  for  analytically 
simulating  parabolic  reflector  antenna  system  before  fabrication, 
installation  and  test. 


The  computer  code,  OSU  PATT,  is  presently  operational  on  the 
RADC-HIS  6180  computer  facility  and  will  be  used  to  analyze  the  performance 
characteristics  of  radar  and  communication  parabolic  reflector  type 
antenna  systems. 


CARMEN  S.  MALAGISI 


Project  Engineer 
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PART  I - THEORY 


CHAPTER  I 
INTRODUCTION 


The  classical  analysis  of  reflector  antennas  was  well  developed 
in  the  late  1940's.  Accordingly,  the  radiation  pattern  in  the  forward 
hemisphere  of  a reflector  antenna  can  be  calculated  by  either  the 
aperture  field  method  or  the  current  distribution  method  [1].  However, 
these  methods  are  rather  slow  for  electrically  large  antennas,  and  in 
general  do  not  predict  the  wide  angle  side  lobes  accurately.  Recently, 
the  wide  angle  side  lobes  were  calculated  [2]  by  the  6TD  [3,4,5]  which 
offers  an  efficient  way  to  obtain  the  radiation  pattern  except  for 
the  forward  axis  region.  Thus  an  overall  pattern  results  from  a 
combination  of  GTD  and  aperture  integration  methods.  Using  this 
approach,  the  far  field  patterns  of  the  circularly  symmetric  parabolic 
reflector  antenna  with  the  feed  at  the  focus  were  calculated  in  [2]. 

The  effect  of  a rapid  field  variation  at  the  edge  of  the  reflector  and 
the  coupling  between  two  reflector  antennas  were  also  analyzed 
using  the  GTD  approach  [6]. 

In  this  report,  the  same  approach  is  extended  to  the  scattering 
from  the  feed  supports.  Also  the  analysis  has  been  extended  to 
include  calculation  of  the  off-principal  plane  patterns.  The  feed 
horn  blockage  is  also  treated  in  a similar  manner  to  that  in  [2]  by 
replacing  the  feed  structure  by  an  equivalent  circular  or  rectangular 
flat  plate  model  whose  area  approximates  the  cross  section  of  the 
feed  structure.  This  approximation  is  justified  in  the  forward 
direction,  where  the  feed  horn  blockage  is  generally  most  significant. 

The  flat  plate  scattering  is  analyzed  by  the  conventional  physical 
optics  approach  [7]. 

On  the  other  hand,  the  feed  support  scattering  is  treated  in  a 
different  way  in  which  the  scattered  field  of  each  individual  strut 
is  obtained  using  the  concept  of  equivalent  current  line  source  [8,9]. 

The  GTD  is  used  to  determine  the  equivalent  line  sources.  The  total 
effect  of  the  feed  support  scattering  is  the  sum  of  the  scattered 
fields  from  each  individual  strut.  Then  the  scattering  from  the  feed 
horn  and  the  feed  supports  is  simply  added  to  the  other  radiation 
components  of  the  antenna.  Thus  the  total  pattern  of  the  reflector 
antenna  is  composed  of  the  radiation  field  from  the  reflector,  the 
direct  feed  pattern  and  the  scattering  from  the  feed  structure. 

' user-oriented  computer  program  is  described  in  Part  II  of  this 
report.  The  input  to  the  program  consists  of  the  E-  and  H-plane  patterns 
of  the  feed,  the  frequency,  the  dimensions  of  the  reflector,  the  positions 
and  diameters  of  the  feed  supports,  and  the  physical  cross  section  of 
the  feed. 
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In  addition  to  the  reflector  antenna  code  described  in  the  body  of 
this  report,  a second  computer  code  developed  by  W.  D.  Burnside  and 
R.  F.  Marhefka  has  been  made  operational  on  the  RADC  Honeywell  computer 
system.  This  code  referred  to  as  the  flat  plate  program  is  used  to  com- 
pute the  far-zone  scattered  fields  for  antennas  radiating  in  the  near 
zone  of  structures  made  of  flat  plates.  In  its  present  form  this  code 
simulates  structures  such  as  buildings,  or  ships  by  a set  of  finite  flat 
plates  forming  a convex  structure  for  which  the  scattering  from  one  flat 
plate  to  another  is  negligible.  Using  the  present  code,  one  can  treat 
the  structure  by  a single  flat  plate,  a rectangular  box,  a rectangular 
pyramid,  etc.  Also  a separate  ground  plane  can  be  introduced.  This 
additional  effort  was  supported  by  the  Naval  Ocean  Systems  Center, 

San  Diego,  California  under  Contract  NOOl 23-76-C-l 371 . 

The  present  code  is  limited  to  one  structure  which  can  be  simulated 
by  as  many  as  14  plates.  This  is  based  on  the  array  dimensions  in  the 
code  and  is  not  a limitation  of  the  theory.  Each  plate  can  consist  of 
6 corners;  however,  each  corner  must  lie  in  a plane  or  the  computer  code 
will  abort.  The  definition  of  the  plates  is  made  by  first  setting  up  a 
fixed  cartesian  coordinate  system  relative  to  the  structure  under  investi- 
gation. The  plates  are,  then,  defined  by  the  location  of  the  corners. 

The  antenna  location  is,  also,  specified  in  the  same  coordinate  frame. 

One  should  note  that  the  fixed  coordinate  system  should  be  chosen  such 
that  one  can  easily  define  the  structure.  The  program  has  the  flexibility 
to  handle  arbitrary  pattern  cuts  relative  to  this  coordinate  system  as  is 
discussed  later. 

The  antenna  presently  considered  in  the  computer  code  is  simulated 
by  a set  of  electric  or  magnetic  elemental  radiators.  There  is  a maximum 
of  six  such  radiators  which  is  limited  by  the  computer  code  dimension 
and  not  the  theory.  Each  electric  or  magnetic  radiator  has  cosine 
distribution,  arbitrary  length,  arbitrary  magnitude  and  phase,  and 
arbitrary  orientation.  This  elemental  antenna  is  considered  initially 
but  can  be  easily  modified  in  that  the  code  is  modular  in  construction. 

In  this  case,  the  SOURCE  subroutines  can  be  easily  exchanged  with  another 
antenna  pattern  subroutine. 

The  present  form  of  the  computer  code  is  not  large  in  terms  of 
computer  storage  and  executes  a pattern  in  short  order.  The  storage  is, 
of  course,  dependent  on  the  dimensions  which  might  vary;  however,  the 
present  code  reguires  approximately  100  K bytes.  It  will  run  a pattern 
cut  of  360  points  for  a flat  plate  structure  with  one  antenna  in  approxi- 
mately 10  seconds  on  a CDC-6600  computer. 
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The  limitations  associated  with  the  computer  code  results  from 
the  basic  nature  of  the  analysis.  The  solution  is  derived  using  the 
Geometrical  Theory  of  Diffraction  (GTD)  technique  which  is  a high  fre- 
quency approach.  In  terms  of  the  scattering  from  a finite  flat  plate, 
this  means  that  each  plate  should  have  edges  at  least  a wavelength 
long.  In  addition,  antenna  elements  should  not  get  closer  than  about  a 
wavelength  to  any  edge.  In  some  cases,  the  previous  wavelength  limit 
can  be  reduced  to  a quarter  wavelength. 

A part  of  the  work  carried  out  under  this  contract  concerned  the 
development  of  diffraction  coefficients  for  perfectly-conducting 
cylindrical  scatterers.  These  diffraction  coefficients  were  used  to 
calculate  the  scattering  from  the  feed  support  of  a reflector  antenna, 
as  described  in  the  body  of  this  report.  However,  they  also  were  used 
to  predict  the  degradation  of  the  pattern  of  the  LAMPS  antenna  caused  by 
the  presence  of  a nearby  cylinder.  The  LAMPS  antenna  is  a 34"  parabolic 
reflector  antenna  with  a nominal  frequency  of  4.6  GHz.  A report  [13]  was 
prepared  on  this  task,  and  since  the  results  were  of  interest  in  the 
analysis  of  a shipboard  antenna  configuration,  its  publication  was 
supported  by  NOSC  through  the  Naval  Regional  Procurement  Office,  Long 
Beach,  California  under  Contract  N00123-76-C-1371 . 
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CHAPTER  II 
RADIATION  PATTERN 


A.  Primary  Feed 

Reflector  antennas  are  classified  according  to  the  geometry  of 
the  reflector  surface  as  well  as  the  shape  of  the  reflector  rim.  In 
this  report,  a focus-fed  paraboloid  with  circular  rim  is  considered. 
However,  this  approach  can  also  be  applied  to  other  kinds  of  re- 
flector antennas.  Since  the  feed  is  located  at  the  focus,  it  is  con- 
venient to  introduce  a spherical  coordinate  system  to  describe  the 
field  of  the  primary  feed  with  origin  at  the  focus  (Pf)  and  the 
y axis  as  the  polar  axis  as  shown  in  Figure  la.  To  describe  the 
radiation  from  the  feed  at  the  reflector,  another  spherical  coordinate 
system  also  with  origin  at  the  focus,  but  having  the  z axis  as  the 
polar  axis  is  used  (see  Figure  Ib^ 
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Figure  1.  Coordinate  systems  for  the  primary  feed. 
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The  relations  between  the  coordinate  variables  and  unit  vectors 
of  these  two  systems  are  given  by 


cos 


i//  = cos  <|>^  sin 


0) 


- cos  (fi'  sin  ij;  = sin  sin 


(2) 


: _ „ sin  41'  cos  1|;  ^ : /COS  41'  ^ 

- 6f  ^ ^ 


- s cos  4.' 

' - - ®f  ~T^ 


‘ sin  41  cos  Ip 
’f ^ 


(3) 

(4) 


where 


2 2 

= Nl-sin  4^sin  4>' 


(5) 


The  coordinate  systems  for  the  reflector  geometry  are  shown  in 
Figure  2 where  (R,e,4>)  are  the  spherical  coordinates  of  the  far  field 


(c)  (d) 

Figure  2.  Coordinate  systems  for  the  reflector. 


point.  If  we  assume  that  the  reflector  is  in  the  far  zone  region  of 
the  feed  source,  the  field  of  the  primary  feed  having  the  same  polari- 
zation in  the  far  zone  as  a dipole  oriented  in  the  y direction  is 
given  by 
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f - -jkR' 

F = - e^A  g(o^.<^^)^-^T — 

, ■ F (5  £.°^  Y^"  ♦.'  - ; Hf^)  f(».«')  (6) 


where  g(6f.<(>f)  = fCi/j,*')  is  the  primary  feed  pattern, 

A is  set  equal  to  F,  the  focal  length  of  the  reflection, 
for  convenience 

2F 

and  R'  = y v co’s'ip  distance  from  the  feed  to  the 

reflector  surface. 

TI 

In  the  off-principal  planes  ((|)  ^ 0,  7),  the  feed  patterns  are  approxi- 
mated by  interpolating  between  E-  ana  ll-plane  patterns  in  the  following 
way. 


f(l^,'i'' ) = ■fj  ■ ■^5  '^os  2(p' 

where 

Tj  - 2 


and 


(7) 


(8) 


fg(l')  = f(l',4>'  = j) 

fh(4-)  = = 0)  . (9) 


B.  Aperture-Integration 

In  this  and  the  next  two  sections  we  are  considering  the 
antenna  pattern  without  aperture  blockage  or  "to  calculate  most  of 
the  scattering  from  the  reflector".  The  6TD  together  w’ uii  the  direct 
feed  radiation  provides  a very  efficient  method  to  calculate  most  of 
the  antenna  pattern.  The  GTD  fields  can  be  calculated  from  only  the 
feed  illumination  of  the  aperture  edges.  However,  the  main  beam  and 
the  first  few  sidelobes  depend  on  the  fields  over  the  entire  aperture. 
Consequently,  the  classical  technique  of  aperture  integration  is  used 
to  calculate  the  pattern  in  the  forward  axial  region. 
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Using  the  concepts  of  geometrical  optics  and  conservation  of 
power,  it  can  be  shown  [2]  that  the  magnitude  of  the  aperture  field 
is  related  to  that  of  the  incident  field  by 


It  follows  from  the  above  expression,  adding  the  appropriate  polari- 
zation and  a phase  factor,  that  the  aperture  field  is  given  by 

-jkRo 

r^(p‘,<^')  = e^  F f(KV)  ^-^T—  (11) 

where 

^r  ” “ F sin4>'  + cos^i'] 

1 1 

= + ^ [j  sin24i'  (1  - cosi|/)  x 

- {cosiii  sirt^ip'  + cos^'p')  y]  . (12) 

Ro  is  the  distance  from  the  feed  to  the  rim  of  the  reflector  (see 
Figure  4) , and 

p ' = R'  sini/i  . (13) 

The  coordinate  system  for  describing  the  aperture  field  is  shown 
in  Figure  2a  in  which  the  origin  is  at  the  center  of  the  reflector 
aperture. 

In  evaluating  the  radiation  from  the  aperture,  a spherical 
coordinate  system,  shown  in  Figure  2b,  and  also  centered  at  0,  is 
introduced  to  describe  the  far  zone  field  of  the  reflector  antenna. 
By  the  eguivalence  principle  with  image  theory,  the  eguivalent 
magnetic  current  induced  by  the  aperture  field  is 

K = 2 X z 

The  pattern  function  of  this  eguivalent  current  source  is  given  by 

F . . ^ ti<(p'.r)  « f']  (14) 

0 0 
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where  r'  is  the  unit  vector  along  r*  , the  distance  between  the  field 
point  and  the  source,  and  a is  the  radius  of  the  aperture. 

In  substituting  ^ into  Equation  (14),  we  consider  the  x and  y 
polarizations  separately.  As  in  [2],  the  aperture  field  E^  can  be 
approximated  by  interpolating  between  the  aperture  distributions  along 
the  principal  axes. 


E^  = Ej(p',f)  = - ; F fj^)  ^ 


E^  = eJ(p'.O)  = - ; F f^(.j;)  ^ 
Thus  the  aperture  field  is  given  by 


Ey(p',<^')  = T^(p')  - 6^(p')  COS2V 

where 


T^(p')  = 


E^  + E^ 
e h 


e"  - 

e h 


6f(p')  = 


Hence  for  the  y component  of  the  aperture  field  in  Equation  (14),  we 
have 


Ky  X r'  = 2(y  E^  x z)  x r' 


2 X E“x(x  sine  cosiji  + y sine  sin<^  + z cose) 


= 2(z  sine  siniji  - y cose)  E‘ 


where  we  have  assumed  r for  far  field  calculation. 
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The  pattern  function  of  Ey  is  therefore  given  by 


Fy(e,<^)  = (y  cose  - z sine  sin(ji)  ^ 


(19) 


where 


2tt  a 


'rll  ^a  gjkp'sinecos(ifi-ii>' ) 


0 0 


The  ())'  integration  can  be  analytically  evaluated  using  the  well-known 
properties  of  the  Bessel  function  Jp;  thus  we  get 


a 

ly  = 27t  J [T^(p')  J^(kp'sine)  + 6^(p')  cos2(kJ2(kp'sine)]p'dp' 
0 

(20) 

The  above  integration  is  numerically  evaluated  by  Simpson's  rule. 

Next  consider  the  x polarized  aperture  field.  From  Equations 
(11)  and  (12) 


E®  = ^ sin2(<>' (1-cos  tk)  F f(i(/,<k')  — 


-jkR, 


(21) 


In  the  principal  planes  (4>  = 0,  7),  these  x polarized  fields  vanish. 
Thus  only  the  y-component  of  the  aperture  field  contributes  to  the 
radiation  patterns  in  the  principal  planes.  In  the  off-principal 
planes,  the  x polarized  component  e9  may  have  significant  contribution 
to  the  radiation  pattern  which  reaches  a maximum  in  the  45°  plane. 

For  the  E^  component  in  Equation  (14)  we  have 


K^x  r'  2(-x  cose  + z sine  C0S4>)  E^  . 


(22) 


Thus  the  pattern  function  Fj(  becomes 

^^ (&.<>')  = (x  cose  - z sine  cos<|i)  ^ 1^ 


(23) 


where 


I 


(24) 


2v  a 
0 ^ 


Using  the  interpolation  for  the  primary  feed  pattern  as  given  by 
Equation  (8),  the  above  integral  becomes 


F e 


-jkR  2ir  a 


j j [sif^2(j>'(l-cosi»>)][fy-fgCos2(j>'] 


X e 


0 0 

jkp-sinecos(<^-<^')p.^p.^^. 


-jkR 

= ttF  e 


0 


2<t»(l-COStj;) 


J2(kp 'sine) 


+ 1 sin4((i(l-cosij;)  ^ J^(kp'sin0)]p'dp'  . (25) 

In  carrying  out  the  ♦ integration,  we  have  made  a stationary  phase 
approximation  that  (<>'  =■  (j>  in  the  expression  for  B,  i.e., 

B = ^1  - sin^ij/  sin2(().  Again  the  above  integration  is  evaluated  by 
Simpson's  rule.  Flence  the  total  far  field  as  obtained  from  aperture 
integration  is  given  by 


p-JkR 

E(R)  = [F^(e,.^)  + Fy(e.<^)]  (26) 

which  is  expressed  in  spherical  components  by 
F"^  = - (e  siniji  + 0 cose  cos(f)  ^ 1^ 
and 

= (e  cos4i  - (|>  cose  sint)  ^ 
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c. 


Two  Point  Method  Using  GTD 


The  wide  angle  side  lobes  of  a circular  reflector  have  previously 
been  analyzed  by  using  the  two  point  method  based  on  GTD  [2].  This 
method  states  that  the  diffracted  field  from  the  paraboloid  is  con- 
tributed mainly  by  two  stationary  points,  Q]  and  Q2  on  the  rim  of  the 
reflector  (see  Figure  3).  The  detailed  formulations  of  the  diffracted 
field  in  the  E-  and  H-planes  are  given  in  [2],  The  off-principal 
plane  diffracted  field  pattern  can  be  obtained  by  using  a similar 
procedure.  The  far  zone  diffracted  fields  from  and  are  given  by 


r^(p)  = E^(Q^)  • D(q^)  I 


a 

sine 


g-jk(R-asine) 

_ 


and 


(27) 


^(P)  = E^Q2) 


°^^2^  ^sTne  R 


(28) 


respectively,  where  D is  the  dyadic  diffraction  coefficient  for  a 
curved  edge,  Et(q-|^2)  is  the  electric  field  of  the  feed  at  Qi  2 
and  a is  the  radius  of  the  reflector  aperture.  For  rays  normally 
incident  on  the  edge,  as  is  the  case  here,  the  dyadic  diffraction 
coefficient  [4]  can  be  expressed  as: 

D = e e + p^  p (29) 


where  Dg  and  are  the  scalar  diffraction  coefficients  for  the  soft 
and  hard  boundary  conditions,  respectively,  and  are  given  in  [4], 
and 

e is  the  unit  vector  tangent  to  the  edge, 

p = e X I;  I being  the  unit  vector  in  the  direction  of  the 
incident  ray, 

P(j  = e X d;  d being  the  unit  vector  in  the  direction  of  the 
diffracted  ray. 

For  the  ray  diffracted  at  Qi 

e = h 

P(j  = e X d = e, 

p = exl^=^xl^ 

til  = o is  the  half  angle  spanned  from 

the  focus  to  the  rim  of  the  reflector. 
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A 


(a)  FRONT  VIEW 


(b)  SIDE  VIEW 


Figure  3.  Geometry  of  the  edge  diffracted  rays 
in  an  arbitrary  <fi  plane. 
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Figure  4.  Different  regions  illuminated  by 
the  edge  diffracted  rays. 

and 

R'  = R is  the  distance  from  the  focus  to 
0 

the  rim  (see  Figures  3 and  4). 

Combining  Equations  (6),  (27)  and  (29)  with  and  (i.e., 
B = Bq  = ^1  - sin^a  sin2()i),  the  diffracted  field  from  point  Qi 
(see  Figure  3)  is  given  by 
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E^(P) 


rf 


E («,♦)  U ^ D^(Q^)  + e (♦  X I^)  D^(Q^)] 
-jk(R-asin0) 


a e 


me  R 


0 

■jkR, 


Rq  ^sine 


a e 


-jk(R-asine) 


Similarly,  the  diffracted  field  from  Q2  is  given  by 

^(P)  = [♦  ^ D^(Q2)  - e D^(Q2)]  f f(a,4.) 


B 

0 

-jkR 


y ® 

""  T 


-1; 


ne 


-jk(R+asine )+j  j 


(30) 


(31) 


The  reflector  rim  is  illuminated  by  the  feed  pattern  as  interpolated 
from  the  E-  and  H-planes.  Thus 


T(ci,(J>)  = f-j-(a)  - fg(o)  cos2<> 


(32) 


where  fj  and  f,5  are  defined  in  Equations  (3a)  and  (3b).  The  diffracted 
field  in  the  different  GTD  regions  is  given  by 


= 


ft  < ■^ 

e 

I t 

®t  - ® 


(33) 


where  0^  is  the  angle  of  the  shadow  boundary  for  Q2  (see  Figure  4). 
Then  the  field  of  the  primary  feed  is  superimposed  on  the  diffracted 
field  in  the  region  from  0=0  to  e = n-a  which  is  the  shadow 
boundary  of  the  field  from  the  feed. 
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D.  Ring  Current  Method 


In  the  rear  axis  region,  the  ring  current  method  is  used  to 
calculate  the  diffracted  field,  since  the  two  point  method  is  not 
valid  there  [2].  The  E-plane  and  H-plane  field  patterns  for  0 close 
to  V have  been  formulated  in  [2]  as 


-j(kR  - J) 

^e  " ■ ® ^ ^ R tAiJo(x)  + A2J2(x)  + A3J^(x)],  (34) 

where  x = kasine, 

^1  “ - Dg  cos  e)  + ^ (D^  + Dj  cos  e)  , 

^2  ~ ~ ^f^^h  * ^s  “ *^f  ^^h  " ^s  ’ 

'^f 

A3  = (D^  + Dj  cos  6)  , (35) 

and 

-j(kR  - J) 

= y ^ ^ R [B^  Jq(x)  + B2  J2(x)  + B3  J4(x)], 

(36) 

where 

^f 

= T^(D^  - cos  0)  - ~ (D^  + cos  0), 

^2  " "^f^^s 

B3  = - ^ (D^  + cos  0).  (37) 

Note  that  the  diffraction  coefficients  Ds  and  are  assumed  to  have 
their  rear  axis  values,  i.e.,  for  0 = it. 

In  the  off-principal  planes,  the  field  in  the  rear  axis  region 
is  obtained  by  interpolating  between  Ee  and  Eh  as  follows: 

Ey  = Eg  sin^<>>  + E^  cos^(|».  (38) 


lb 


This  can  be  transformed  to  spherical  components  by 


Eg  = Sin<t> 

E^  = Ey  coS(ji  . 

since  ir-e  is  assumed  to  be  small. 


(39) 

(40) 
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CHAPTER  III 

FEED  STRUCTURE  SCATTERING 


In  practice,  the  scatterinq  from  the  feed  horn  and  the  feed 
supports  increases  the  side  lobe  levels  and  reduces  the  gain  of  the 
antenna.  Therefore,  the  blocking  effect  of  the  feed  structure  should 
be  taken  into  account. 

To  analyze  the  scattering  from  the  feed  structure,  two  different 
models  are  used.  An  equivalent  line  source  model  is  used  for  feed 
support  scattering  and  a flat  plate  model  is  used  for  feed  horn 
scattering. 

A.  Equivalent  Line  Source  Model 
F'or  Teed'  Support  Scattering' 

The  effect  of  the  feed  support  blocking  is  usually  estimated 
by  the  projected  shadow  method  in  which  the  feed  support  or  strut  is 
replaced  by  its  effective  shadow  on  the  reflector  plane.  The 
radiation  field  from  the  shadow  area  is  then  calculated  and  subtracted 
from  the  antenna  pattern  [10,11].  In  this  report,  an  alternative 
approach  is  developed  in  which  the  scattered  fields  from  the  struts 
are  computed  by  the  equivalent  current  approach  [7].  Since  struts 
often  take  the  form  of  circular  cylinders,  only  metallic  circular 
cylinder  struts  are  considered  here.  Also,  the  incident  field  on  the 
struts  is  assumed  to  be  the  reflected  wave  from  the  reflector  only. 

That  is,  the  interaction  of  the  direct  feed  radiation  and  the  strut  is 
not  considered.  The  strut  scattered  field  as  reflected  by  the  re- 
flector is  also  neglected.  Therefore,  the  feed  support  scattering 
problem  is  equivalent  to  a circular  cylinder  scatterinq  with  a plane 
wave  incident  at  a specified  incident  angle;  however  the  amplitude 
of  the  plane  wave  can  be  a slowly-varying  function  of  position.  The 
equivalent  current  approach  is  used  to  take  into  account  the  vari- 
ations of  the  incident  field  along  each  strut,  which  includes  the 
effect  of  the  finite  length  of  each  strut.  In  this  approach,  the 
scattering  from  each  element  of  a strut  is  assumed  to  be  the  same  as 
that  from  an  element  on  an  infinite  cylinder  with  the  same  incident 
field. 

Let  (x,n,C),  (c,Y,4)  and  (r,a,y)  denote  the  rectangular,  cyindrical 
and  spherical  coordinates  of  the  strut  system,  respectively.  Consider  a plane 
wave  incident  upon  an  infinite  conducting  cylinder  of  radius  a,  located  at 
the  origin,  with  an  incident  angle  8,  at  y'  = tt,  as  shown  in  Figure  5. 
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X 


Figure  5.  Geometry  of  the  scattering  cylinder  with  an  incidence 
angle  e. 

(Note  that  and  both  lie  in  the  plane  of  incidence, 
which  contains  the  incident  ray  and  the  cylinder  axis.) 


The  incident  field  can  be  expressed  as 

F'  = Q cose  + I sine)  E^  ^ ^ <^°se)  ^ (4^) 

The  field  component  parallel  to  the  axis  of  the  cylinder  is  given  by 
Ej  = E,,  sine  sine  + C cose)  . (42) 

Noting  that 

X = ^ COSY  (43) 


.1 


we  have 


e;  . E„  sinB  e-J'" 


(44) 


For  a wave  incident  at  an  angle  3 with  respect  to  the  cylinder  axis, 
the  analysis  of  the  cylinder  scattering  closely  follows  that  for 
3 = 90°  [7].  Thus 

Ej  = E„  sin3  ^ I j""  J^(kc  sin3)  . 


Then  the  scattered  field  component  from  the  infinite  circular  cylinder 
is  given  by 

E^  = E„  Sins  ^ f j‘"  H^^^kcsinS)  x e^"^  . 


In  the  far  field  zone 


n n 


(45) 


E^^  » E„  sins 


-~i-  ~ I a e 

TflTcsinS  n 


jny 


Em  sins 


j-jks 


sinS 


I I a„  cos  ny 
nk  n n 


where 


J (ka  sin  ) 

^n  = - 


'(ka  sin  ) 
n ' 


^n 


n = 0 
otherwi se 


and 


s = t/sinS 


The  scattered  field  for  parallel  polarization  can  be  represented  by 
E^  -jks 
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where  the  diffraction  coefficient  is  given  by 


D^(Y.e) 


I a„  cosny 
n=0 


(47) 


Comparing  with  the  field  of  an  infinite  line  source  with  current 
I = g-jkCcos3^  which  is  given  by 


E = 


a 


k^I 

0 rs 

Awe  iliik 


(48) 


we  obtain  an  equivalent  electric  current  line  source  for  the  con- 
ducting cylinder  with  a current  equal  to 


|¥ 


(49) 


Then  the  scattering  from  a finite  section  of  conducting  cylinder  can 
be  obtained  by  calculating  the  radiation  from  the  equivalent  finite 
length  line  source.  Thus  the  equivalent  current  I for  the  parallel 
polarized  incident  field  becomes 

1(5)  « I„(5)  (50) 

and  the  far  zone  radiation  field  is  given  by 


E 


a 


D5(Y.e) 


T||(C) 


gjkC (cosa+cos3) 


dC 

(51) 


where  (r.a.y)  specifies  the  spherical  coordinates  of  the  strut  system. 
Similarly  for  the  electric  field  component  perpendicular  to  the  axis 
of  the  cylinder,  the  scattered  field  can  be  obtained  from  H„  as 


H 


a 


Dh(’^*6) 


r 


sina 


H„(£) 


^jk? (COSO+COS0) 


d^ 

(52) 
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where 


D 


h 


(y,0) 


I En  b cosny 
n=0  " 


and 


J^(ka  sine) 

" sine) 

The  corresponding  electric  field  for  the  perpendicularly  polarized  incident 
field  is  given  by 


E = 
Y 


D^(y>0) 


sina 


X 


Ej_{0  de 


(53) 


B.  Coordinate  Transformation 


In  order  to  add  the  fields  scattered  from  the  feed  struts  to 
the  radiation  field  from  the  reflector,  each  strut  coordinate  system 
needs  to  be  transformed  into  the  reflector  coordinate  system.  As 
described  in  the  previous  section,  the  rectangular  coordinates  of  the 
strut  coordinate  system  are  denoted  bv  (x.n,4);  and  (r,a,y)  denotes 
the  spherical  coordinates.  The  relationship  between  their 
respective  unit  vectors  is  given  by 


“ ^ - 

r 

X 

a 

= [A] 

n 

f— 

> 

L 

r 

1 

where 


[A] 


sina  cosy 
COSa  cosy 
-siny 


sina  siny 
cosa  siny 
cosy 


cosa 

-sina 

0 


(54) 


(55) 


is  the  transformation  matrix. 
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Similarly  for  the  reflector  coordinate  system 


R 

0 

= [B] 

X 

y 

1 

■©-  > 
1 

z 

where 


sine  cos<ji 

si  ne 

si  n(ji 

cose 

[B]  = 

cose  cos4i 

cose 

sin(}> 

-sine 

-si  n$ 

COS<() 

0 

(57) 


Now  let  us  transform  the  (x,y,z)  coordinates  into  the  (x.n.C)  ones. 
The  transfonridtion  consists  of  two  steps;  first  rotate  the  x axis  an 
angle  about  tho  z axis,  then  the  z axis  is  rotated  an  angle  Sq 
about  the  y'  (or  n)  axis,  as  shown  in  Figure  6.  The  resultant  trans- 
formation can  be  expressed  by  the  following  representation: 


cose„  cosiii„ 
0 0 

COS0^ 

-sine 

-sin^Q 

COS<fro 

0 

sine^  '^os^Iq 

sine  Sint 
0 0 

cos6o 

(58) 


If  we  denote  the  angle  between  the  strut  along  the  t axis  and  the 
negative  z-axis  as  e,  and  the  angle  between  the  x-axis  and  the  strut 
projection  on  the  xy  plane  as  (see  Figure  7),  we  have 

% = (59) 


(60) 


and 


1 

> X > 

J 

= [T] 

X 

y 

LiJ 

z 

where 


[T] 


r~ 


-cose  COSi}>^ 

-cose  siniji^ 

-sine 

-sin4>^ 

COS(li^ 

0 

sing  cosil>^ 

sine  sini^^ 

-cose 

(61) 


(62) 


(a)  ROTATE  4>o  ANGLE  ABOUT  Z AXIS 


Figure  6.  Coordinate  transformations  between  the  strut 
and  the  reflector  coordinate  systems. 


Then  the  relation  between  the  two  spherical  coordinates  (R,e,4i)  and 
(r.a.y)  is  given  by 


< 

L_ 

R 

^ > c > 

1 

= [A][T][B]-1 

e 

(63) 


where  [B]~^  is  the  inverse  matrix  of  [B],  To  express  a and  y in 
terms  of  and  415,  the  components  of  x,n,5  are  written  in  the 


(x,y,z)  coordinate  system  as 

X = - X cose  cosi)>^  - y cose  sim^^  - z sine  (64) 
n = - X sini}>^  + y cos(Ji^  (65) 
i = X sine  cos4i^  + y sine  sin<^^  - z cose  . (66) 

But 

X = R sina  cosy  (67a) 
n = R sina  siny  (67b) 
5 = R cosa  (67c) 

and 

X = R sine  cos<()  (68a) 
y = R sine  sin()>  (68b) 
z = R cose  (68c) 


Putting  Equations  (67a)  and  (68)  into  Equation  (64)  and  combining 
terms,  we  obtain 

sina  cosy  = - [sine  cose  + cose  sine]  . (69) 

Similarly,  putting  Equations  (67b)  and  (68)  into  Equation  (65),  we  get 
sina  siny  = sine  sin (1(1-41^)  . (70) 

Dividing  Equation  (70)  by  Equation  (69),  we  have 
sine  sin(4>-4>j) 

^ “ -[sine  cose  cos(4'-4'5)  + cose  sineT  ' 
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Also,  substituting  Equations  (67c)  and  (68)  into  Equation  (66),  we 
get 


cosa  = sine  sin0  cos{<^-(p^)  - cose  cose  . 

(72) 

For  a given  observation  direction  e and  the  angular  coordinate 

a and  y for  a given  strut  can  be  calculated  from  Equations  (71)  and 
(72),  where  6 and  specify  the  orientation  of  the  strut.  The 

field  components  in  the  strut  coordinate  system  can  be  transformed 
to  the  reflector  coordinate  system  by  using  the  following  dot 
products  obtained  from  Equation  (63)  as 

a-e  = - cose  cos(ifi-i)>^)[sina  sine  + cosa  cosy 

cose] 

+ COSO  siny  sin(<j)-<ji^)  cos0 
+ sin0(cosa  cosy  sine  -sina  cose) 

(73a) 

a*(|i  = [sina  sine^  + cosa  cosy  cose]  sin(4i-4i^) 

+ cosa  sin  y cos  (i)>-()>^) 

(73b) 

y*e  = siny[cose  cose  cos(4>-4ig)  - sine  sine] 

+ cosy  cose  sin((j)-<{i^) 

(73c) 

= - siny  cose  sin(4>-(fi^)  + cosy  cos(iJ>-4i^) 

(73d) 

Thus  the  far  field  scattered  from  the  strut  can  be  expressed  in  the 
(R,e,i)>)  coordinate  system  as 

E„=e*aE  +e*yE 
e a y 

(74) 

E.  = ^-a  E + ()i*y  E 
<t>  0 ^ y 

(75) 

E^  and  E^  are  given  by  Equation  (51)  and  (53),  respectively  in  which 
E„  = eJ  sin0^  (76) 

and 

Ei.  = Ey  cosiji^  (77) 

as  shown  in  Figure  7. 
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In  addition  to  the  coordinate  rotation  described  above,  the 
transformation  also  consists  of  a linear  displacement  of  the 
coordinate  centers  along  the  z-axis  as  shown  in  Figure  8.  Thus 


Figure  8.  Linear  axis  placement  of  the  strut  and 
the  reflector  coordinate  systems. 


r=R-ZQCOs6  in  the  phase  expression  for  the  far  field  approximation 
in  Equations  (51)  and  (53),  and  r=R  in  magnitude.  The  quantity  Zq 
is  the  distance  between  the  two  coordinate  centers,  i.e.,  the  distance 
from  the  center  of  the  reflector  aperture  to  the  intersection  of  the 
strut  on  the  z-axis. 


C.  Feed  Horn  Scattering  Model 


A conmonly  used  method  to  analyze  the  feed  horn  blockage  is  to 
assume  that  the  scattering  from  the  feed  aperture  is  the  same  as 
that  from  a conducting  flat  plate  whose  area  is  equal  to  the  cross 
section  of  the  feed.  This  approximation  is  justified  in  the  forward 
region,  where  the  feed  blockage  is  generally  most  significant.  Then 
the  effect  of  the  aperture  radiation  being  blocked  by  the  feed  can  be 
approximated  by  adding  the  pattern  of  the  equivalent  flat  plate  to 
that  of  the  reflector.  Since  the  feed  aperture  is  relatively  small 
compared  with  the  antenna  aperture  and  is  located  at  the  center,  the 
aperture  field  can  be  assumed  uniform  over  the  flat  plate.  So 


(78) 
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(b)  RECTANGULAR 


Figure  9.  Flat  plate  models  for  the  feed  horn  scattering. 

for  a feed  with  circular  aperture  of  radius  C (see  Figure  9a),  such 
as  a circular  feed  horn.  The  far  zone  scattered  field  for  small  6 is 
given  by 


-jkR 
Tl 


2nc 


J-|  (k  c sine) 
k c sine 


(79) 


If  the  feed  is  a rectangular  waveguide  or  horn,  the  scattered  field 
is  given  by 
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y J 


. e 


-jkR 


XR 


sin(|^  sine  sin(<>) 
^ sine  sin(|i 


ab 


sin(^  sine  cos^) 


~ sine  cos4>) 


(80) 


where  a and  b are  the  width  and  height  of  the  flat  plate  model, 
respectively,  {see  Figure  9b).  The  scattered  field  components  can 
be  obtained  in  spherical  coordinates  by  using  Equations  (39)  and  (40). 
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CHAPTER  IV 

RESULTS  AND  DISCUSSION 


Previous  analyses  [10,11,12]  of  feed  support  scattering  have 
usually  been  limited  to  angles  near  the  main  beam  region.  However, 
the  cylinder  scattering  approach  developed  in  this  report  provides  a 
model  to  treat  the  wide  angle  scattering  of  each  strut.  Therefore, 
in  addition  to  the  gain  loss,  this  approach  also  gives  the  effect  of 
strut  scattering  on  the  wide  angle  side  lobes  of  the  reflector. 

According  to  the  geometrical  optics  assumption,  the  forward 
scattering  from  a cylinder  is  identical  to  that  from  a strip  with  the 
same  cross  section  provided  that  the  diameter  of  the  cylinder  is 
large  enough.  Thus  the  field  obtained  from  the  two  solutions  should 
be  nearly  the  same  in  the  forward  region.  In  Figure  10,  the  scattered 
pattern  of  a rectangular  strip  is  compared  with  that  of  a vertical 
circular  cylinder  (g  = 90°)  with  uniform  incident  field,  whose  E-field 
is  parallel  to  the  cylinder  axis.  As  shown  by  the  curves,  the  field 
near  the  forward  axis  in  the  H-plane  (transverse  to  the  cylinder  axis) 
agrees  within  1 dB  for  kap  = 7,  where  ap  is  the  radius  of  the  strut. 
For  larger  kap,  the  difference  is  even  smaller  since  the  geometrical 
optics  approximation  of  the  rectangular  strip  for  the  cylinder  gets 
better.  For  wide  angles  (e  > 20°),  the  agreement  becomes  poorer, 
because  the  shadow  aperture  model  is  not  valid  there.  Another  com- 
parison is  made  between  the  Induced  Field  Ratio  (IFR)  coefficients  in 
the  paper  by  Rusch  et  al  [12],  and  those  calculated  from  the  equations 
in  this  report.  The  IFR  coefficient  is  defined  to  be  the  ratio  of  the 
field  forward  scattered  by  the  cylinder  to  that  of  a flat  strip  with 
the  same  physical  cross  section  as  the  cylinder.  Thus  the  IFR  can  be 
obtained  from  the  ratio  of  Equation  (51)  with  E||  = 1 , = 0 and 

^2  = b to  Equation  (80)  with  a = 2ap  and  b = L.  Values  of  the  IFR 
obtained  from  this  ratio  are  compared  in  Table  1 with  those  calculated 
by  Rusch  et  al , where  the  IFR  coefficients  of  [12]  are  taken  from 
Figure  2 of  that  paper.  The  consistency  of  the  two  methods  in  the 
forward  axis  region  is  evident.  Since  the  IFR  as  defined  in  [12]  is 
restricted  to  forward  scatter,  the  results  of  [12]  cannot  be  used  in 
this  work,  where  the  bistatic  scattering  from  the  cylindrical  strut 
is  required. 
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STRUT  DIMENSION  FLAT  STRIP  DIMENSION 


Comparison  of  the  scattering  from  a circular  cylinder 
(vertical  strut)  and  a flat  strip  in  the  H-plane. 


Table  1 

IFR  Coefficients  for  a Circular  Cylinder 


Arq  = kasinS 
X = Arq/TT 


IFR  (E-PLANE) 

IFR  (H-PLANE) 

Arg 

B 

Eqs.  (51), (80) 

Ref. [12] 

Eqs. 

(51), (80) 

Ref. [12] 

2.54 

1.34  /-19.8° 

1.33  /-20.3° 

0.78 

/22.9° 

0.79  /21.9° 

4.62 

' 1.47 

1.22  /-14.3° 

1.22  /-15.0° 

0.85 

/15.7° 

0.85  /16.4° 

6.74 

2.15 

1.16  /-11.5° 

1.16  /-12.1° 

0.88 

/12.4° 

0.89  /11.5° 

7.02 

2.24 

1.16  /-11.2° 

1.15  /-11.7° 

0.89 

/ll .5° 

0.90  /ll .0° 

Fiqure  11  shows  the  geometi^y  of  a practical  reflector  antenna 
system  which  contains  three  struts  with  diameter  2an  = 0.37",  e = 68° 
and  Iks  = 90°,  210°,  and  330°,  respectively.  The  diameter  of  the 
reflector  aperture  is  24",  with  F/D  = 1/3.  The  rectangular  feed  horn 
is  mounted  in  a circular  flat  plate  with  radius  c = 1.2".  The 
H-plane  pattern  of  this  antenna  system  for  a frequency  = 11  GHz  is 
shown  in  Figure  12.  The  pattern  for  the  principal  total  field 
corresponds  to  the  principal  polarization  and  consists  of  the  total 
radiation,  including  that  from  the  reflector,  the  direct  feed  pattern, 
the  feed  model  scattering  and  the  feed  strut  scattering.  Also 
shown  in  Figure  12  is  the  pattern  for  the  radiation  from  the 
reflector  and  the  direct  feed  radiation,  without  the  feed  system 
scattering. 

As  seen  from  Figure  12,  the  scattering  from  the  feed  and  feed 
struts  causes  little  effect  on  the  main  beam.  This  is  predictable 
since  the  radiation  field  from  the  reflector  itself  dominates  in  this 
region,  as  can  be  seen  in  Figure  13  which  shows  the  feed  scattering 
components  separately.  Beyond  the  main  beam  and  the  first  few  side 
lobes,  the  scattering  from  the  struts  and  the  feed  dominates  and 
substantially  raises  the  side  lobe  levels  for  angles  less  than  35°. 
Figure  14  shows  the  results  for  the  off-principal  plane  <p  = - 15°  in 
which  the  strut  scattering  is  high  for  angles  out  to  80°.  This 
occurs  because  the  scattering  cone  (see  Figure  15)  of  the  strut  at 
Iks  = 90°  has  a maximum  near  0 = 65°  in  the  ik  = - 15°  plane,  in 
addition  to  the  maximum  on  the  forward  axis.  The  other  two  struts 
make  less  contribution  to  the  scattered  field  in  this  plane  since 
their  scattering  cones  are  located  in  different  directions.  Therefore 
the  feed  support  scattering  effect  in  the  <k  = - 15°  plane  comes  mainly 
from  the  strut  at  ikj  = 90°.  The  scattered  field  patterns  of  the 
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Figure  12.  H-plane  pattern  (i|i=0)  of  the  reflector  antenna. 
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Figure  13.  H-plane  pattern  (4'=0)  of  the  feed  strut  scattering  and 
feed  horn  scattering. 


Figure  14.  Off-principal  plane  pattern  ((|>=-15°)  of  the  reflector  antenna. 


Figure  15.  Scattering  cone  of  a cylinder. 


struts  and  the  feed  horn  are  shown  in  Figure  16.  Note  that  the  total 
pattern  is  no  longer  symmetric  with  respect  to  the  z-axis  for  most 
off-principal  plane  cuts  as  can  be  seen  by  comparing  Figure  14  and 
Figure  17  for  the  (p  = - 15°  and  41  = + 165°  plane  pattern.  Even 
though  the  feed  pattern  (and  thus  the  aperture  field)  is  assumed  to 
be  symmetric,  the  strut  geometry  is  not  symmetric  in  this  plane. 

The  antenna  patterns  in  these  figures  correspond  to  the 
principal  polarization,  defined  in  terms  of  a Huygen's  source.  Thus 
the  6 and  ♦ components  can  be  transformed  to  principal  and  cross 
polarized  components  by 

^rinc  = ^ 

^cross  = ^6  - 

The  fields  scattered  from  the  feed  and  struts  are  shadowed  by  the 
reflector  surface  for  much  of  the  rear  hemisphere.  However,  the 
shadow  boundary  for  the  feed  support  scattering  is  difficult  to 
define  since  part  of  the  scattered  field  is  blocked  by  the  reflector 
as  d approaches  90°,  and  this  blocking  keeps  on  increasing  as  6 passes 
90°  until  the  scattered  field  is  totally  blocked  out  at  e = ti  -a,  the 
shadow  boundary  for  direct  feed  radiation.  For  simplicity,  the 
shadow  boundary  for  the  feed  support  and  feed  horn  scattering  is  set 
at  e = 90°.  This  approximation  should  be  reasonable  since  the  scat- 
tered field  from  the  struts  is  usually  small  at  pattern  angles  6 90°. 

However,  a slight  discontinuity  still  exists  in  the  pattern  at  e = 90° 
because  of  this  truncation. 
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cattering  and  the  feed  horn  scattering. 


Figure  17.  Off-principal  plane  pattern  (4i=+165°)  of  the  reflector  antenna. 


CHAPTER  V 
CONCLUSIONS 


In  addition  to  the  gain  loss,  which  has  been  discussed  in 
previously  reported  analyses  of  feed  strut  scattering  referenced 
in  Chapter  IV,  the  cylinder  scattering  model  developed  in  this  report 
gives  the  effect  of  strut  scattering  on  the  wide  angle  side  lobes 
of  the  reflector  antenna.  The  scattering  from  the  feed  struts  always 
has  its  greatest  effect  near  the  scattering  cones  for  each  strut. 

These  scattering  cones  usually  give  rise  to  a maximum  effect  on  the 
antenna  sidelobes  in  certain  off-principal  plane  patterns.  For 
example,  the  maximum  effect  of  one  of  the  antenna  struts  in  Figure  11 
occurs  in  the  off-principal  plane  with  = - 15°  as  shown  in  Figures 
14  and  16.  When  the  other  contributions  to  the  radiated  field  are 
included,  the  analysis  given  here  can  be  used  to  compute  the  complete 
pattern  of  any  arbitrary  plane  cut  for  a practical  reflector  antenna 
system.  The  basic  radiation  pattern  from  the  reflector  itself  is 
obtained  by  using  GTD  for  most  of  the  pattern  and  aperture  integration 
for  the  main  beam  region. 

Although  only  circularly  symmetric  reflector  systems  with 
focused  feeds  are  discussed  here,  the  approach  can  be  extended  to 
reflector  antennas  with  arbitrary  rim  shapes  and  to  off-set  fed 
systems  as  well.  Once  the  aperture  field  is  known,  the  aperture 
integration  can  be  carried  out  and  the  incident  fields  on  the  struts 
are  available.  Thus  the  scattering  from  each  strut  is  readily 
obtained.  The  same  idea  can  also  be  applied  to  the  struts  of 
rectangular  or  other  cross  sections;  however,  the  diffraction  coef- 
ficients are  more  difficult  to  calculate  in  these  cases. 
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PART  II  - USER'S  MANUAL  FOR  COMPUTER  PROGRAM 
CIRCULAR  REFLECTOR  ANTENNA  WITH  FEED  AND  STRUT  SCATTER"  (OSUPATT) 


A.  Description  of  Computer  Program 

OSUPATT  is  a computer  program  which  calculates  the  total  field 
radiated  by  a circularly  symmetric  parabolic  reflector  including  effects 
of  scattering  by  a number  of  circularly  cylindrical  struts  and  blockage 
by  the  feed  horn.  The  user  is  free  to  specify  the  geometry  of  the 
reflector,  the  size  of  the  feed  horn  cross-section,  both  the  E-  and 
H-planes  of  the  feed  pattern,  and  the  size  and  orientation  of  the 
struts.  This  brief  user's  document  will  indicate  the  methods  used  and 
the  limitations  of  the  program,  and  required  formats  for  input  data; 
details  of  the  theory  are  given  in  Part  I of  this  report. 

The  basic  geometry  used  in  the  program  is  shown  in  Figure  1. 


y 


Figure  1.  Geometry  of  Reflector. 
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The  reflector  points  in  the  +z  direction  with  the  feed  located  at  the 
focus  of  the  parabola.  The  polarization  of  the  feed  is  linear  in  the 
+y  direction;  thus  is  the  principal  polarization  in  the  <^  = 0 plane 
and  Ee  is  the  principal  polarization  in  the  = 90°  plane.  The  E-plane 
of  the  feed  horn  and  the  reflector  is  thus  the  yz  plane  and  the  H-plane 
is  the  xz  plane.  The  struts  must  be  oriented  so  that  each  strut  projects 
through  the  z axis.  At  this  time  the  principal  plane  feed  patterns  must 
be  sytmetric.  These  patterns  are  specified  by  listing  field  values 
for  discrete  values  of  ip,  the  angle  away  from  the  -z  axis  (see  Figure  2). 


Figure  2.  Feed  horn  geometry. 

The  following  sections  describe  specific  aspects  of  the  problem. 


DESCRIPTION  OF  INPUT  PATTERN 


As  mentioned  before,  the  feed  pattern  is  described  in  its 
principal  planes  by  supplying  tabulated  field  values.  The  feed 
pattern  is  approximated  by  a piecewise  linear  function.  A provision 
has  been  added  so  that  the  assumed  field  values  past  the  last  tabulated 
value  will  be  0.  To  be  safe,  one  should  be  sure  to  specify  the  pattern 
explicitly  from  il-  = 0 to  180°. 

Between  principal  planes,  the  field  values  are  simply  interpo- 
lated from  the  values  in  the  principal  planes.  The  E-plane  values  are 
weighted  by  sin^iji  and  the  H-plane  value  by  cos^^i. 
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REFLECTOR  PATTERN  CALCULATION 


The  reflector's  secondary  pattern  is  calculated  by  two  different 
methods  over  two  different  ranges  of  polar  angle  6.  For  small 
kD 

6 (-y  sin  e < 24),  the  aperture  fields  are  integrated  as  specified 

by  physical  optics.  Due  to  the  symetry  of  this  problem,  the  inte- 
gration is  still  one  dimensional  and  thus  relatively  large  apertures 
can  be  used  without  prohibitive  computer  costs.  For  6 angles  past  the 
limit,  the  Geometrical  Theory  of  Diffraction  (GTD)  is  used.  The  fields 
calculated  by  the  program  are  all  normalized  by  the  peak  field  of 
the  reflector  by  itself.  Any  direct  field  from  the  feed  is  also 
included  in  the  total  field. 


FEED  BLOCKAGE 


The  effect  of  blockage  by  a feed  cross-section  is  very  simply 
modeled  in  the  forward  direction  by  Kirchkoff's  method.  When  a 
rectangular  feed  is  specified,  the  scattered  pattern  has  the  familiar 
sin(x)/x  shaped  pattern  in  each  principal  direction  while  for  a 
circular  feed,  the  pattern  has  J](x)/x  variation. 


STRUT  BLOCKAGE 


Scattering  by  circular  cylindrical  struts  is  calculated  by 
first  finding  the  equivalent  currents  on  segments  of  the  struts  and 
then  integrating  the  reradiated  fields  which  are  produced  by  these 
currents.  The  equivalent  currents  are  calculated  from  the  diffraction 
coefficients  of  the  cylindrical  strut  cross-section  and  the  incident 
aperture  ^ield.  The  diffraction  coefficients  are  a summation  of 
eigenmodts  for  the  circular  cylinder. 
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LIMITATIONS 


Due  to  limitations  in  theory  or  the  arrangement  of  the  computer 
program,  the  following  limits  should  be  observed. 

STRUTS 

STRUT  must  not  extend  past  z axis,  i.e.,  RL02  must  not 
RS22 

exceed  . The  case  of  a strut  passing 

through  the  axis  can  be  done  by  using  two  separate 
struts  that  meet  at  the  axis. 

STRUT  diameter  should  not  exceed  four  wavelengths. 

Maximum  number  of  struts  is  4.  Beta  should 
be  in  the  range  0°  to  90°. 

For  best  results,  strut  angle  BETA  should  be 
large  enough  so  that  sampling  points  are  less 

than  a 1/4  wavelength  apart,  i.e.,  sin  3 > . 

FEED  PATTERN 

The  maximum  number  of  points  to  fit  the 
pattern  (N2)  is  50. 


REFLECTOR 


Maximum  diameter  is  100  wavelengths. 


INPUT  FORMAT 


GENERAL  COWENTS 

All  distance  parameters  required  for  input  are  in  inches; 
all  angles  should  be  specified  in  degrees. 

All  decimal  numbers  may  be  anywhere  within  their  specified 
field;  all  integers  must  be  right  adjusted.  The  decimal  point  must 
explicitly  appear  for  all  decimal  numbers. 

Comnents  may  be  typed  on  any  card  to  the  right  of  specified 
fields. 
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REFLECTOR  GEOMETRY  CARD 


PURPOSE:  To  enter  geometry  of  reflector 


PARAMETERS;  (I  integer,  F decimal  number) 

D2  (F)  -diameter  of  reflector 

FREQ  (F)-  frequency  in  MHZ 

FDR  (F)  - ratio  of  focal  distance  f to 
diameter  D 

NAP  (I)  - number  of  aperture  integration  points. 

If  this  parameter  is  less  than  50, 
or  if  it  is  omitted,  50  will  be  used. 
If  over  99  is  used,  99  will  be  used. 


NOTES: 

• The  NAP  parameter  allows  an  increase  of  inteqration 
points  over  both  the  aperture  and  the  strut  for 
large  sized  reflectors,  while  allowing  a smaller 
number  of  points  and  smaller  computer  costs  for 
smaller  problems.  The  value  99  may  always  be  used 
if  the  computer  cost  is  not  objectionable. 
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FEED  PATTERN  CARD 


PURPOSE:  To  provide  data  on  number  of  points  used  to  model 
primary  feed  pattern  . 


CARD: 


PARAMETERS:  (1  integer,  F decimal  number) 

N2  (I)  - number  of  data  points  used  to  determine 
pattern  in  piecewise  linear  function. 


NOTES : 

• This  card  determines  how  many  cards  of  the  next 
format  will  be  read. 
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FEED  PATTERN  DATA  CARDS 

PURPOSE:  To  provide  tabulated  val  ues  of  primary  pattern. 

CARD: 


PARAMETERS:  (I  integer,  F decimal  number) 

PSl  (F)  : Angle  from  axis  of  horn  (-z  axis 

of  reflector,  see  fiqure  2)  where 
tabulated  value  of  field  is  to  be 
supplied . 

FIELD  (F)  : Field  value  of  horn  at  that  PSI  angle. 

NOTES : 

• There  should  be  2 groups  of  cards  of  this  format. 

M2  ^Tabulated  values  to  be  interpolated  by  piece- 
cards  (_wise  linear  function  for  E plane 

N2  TTabulated  values  to  be  interpolated  by  piece- 
cards"^wise  linear  function  for  H plane 

• For  PSI  larger  than  the  last  given  value  of  PSI, 
a zero  is  assumed  for  the  field  value 

• The  maximum  value  of  psi  should  subtend  the  edqe 
of  the  reflector  or  erroneous  results  will  occur. 
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NUMBER  OF  STRUT  PROBLEMS  CARD 


PURPOSE:  To  indicate  to  the  program  how  many  problems  will 

be  done  with  the  same  reflector  and  feed  pattern 
(above  cards)  but  with  different  blockaae  para- 
meters (below  cards) 


CARD: 

A 

( NCASE 


Th«-  nunbors  along  the  top  refer  to  toe  lost  colur.n  in  each  field. 


PARAMETERS:  (I  integer,  F decimal  number) 

NCASE  (I)  - number  of  different  strut  problems 
to  be  done  with  same  reflector  and 
feed  pattern 


NOTES: 

• If  NCASE  is  not  1,  the  entire  group  of  cards  to 
be  described  below  must  be  repeated  NCASE  times 

• If  the  NCASE  field  is  left  blanlc,  NCASE  will  be 
assigned  1.  However,  this  card  must  always  be 
present  even  if  it  is  blan)<. 
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FEED  SCATTERING  MODEL  CARD 


PURPOSE:  To  describe  size  of  feed  aperture  for  blockage 

purposes . 


CARD: 


PARAMETERS:  (I  integer,  F decimal  number) 

IG  (I)  - Code  variable  that  indicates  shape 
of  feed  aperture. 

The  values  are:  1 - feed  is  circular 

2 - feed  is  rectangular 

DIM  (1)  (F)  - These  parameters  are  dimensions 

DIM  (2)  (F)  of  feed  (see  figure  3) 

If  KG  = 1,  DIM  (1)  is  radius 
DIM  (2)  is  ignored 

If  KG  = 2,  DIM  (1)  is  length  of  feed 
in  X direction 
DIM  (2)  is  length  of  feed 
in  y direction 
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COMMON  STRUT  DATA  CARD 


PURPOSE:  To  provide  data  that  is  conunon  to  all  struts. 


CARD: 


/ 

20 

DP2 

BETA 

NS 

The 

nuo)b<?rs  a 

ong 

Che  cop  refer  co  the  Ijsc  colur.n  in  each  field. 

PARAMETERS:  (I  integer,  F decimal  number) 

DP2  (F)  - diameter  of  struts 

BETA  (F)  - angle  of  struts  from  -z  axis 
(see  figure  4) 

NS  (I)  - number  of  struts  (cannot  be  more 
than  4) 


NOTES : 


This  card  causes  NS  cards  of  the  next  format 
to  be  read. 
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INDIVIDUAL  STRUT  DATA  CARDS 


PURPOSE:  To  describe  data  that  may  be  different  for  each 

strut . 


CARD: 


PARAMETERS:  (I  integer,  F decimal  number) 

RS22  (F)  - Radial  distance  of  one  end  of  the 
strut  from  z axis  (see  figure  4) 

ZS22  (F)  - z coordinate  of  this  same  end 

RL02  (F)  - length  of  strut 

PHIS2  (F)  - 4>(phi)  angle  of  projection  of  strut 
in  the  xy  plane 

NOTES : 

• For  good  accuracy  there  is  a relationship 
between  NAP,  diameter  of  reflector  in  wave 
lengths,  and  angle  BETA.  For  a given  NAP, 
and  diameter,  the  more  nearly  parallel  the 
struts  are  to  the  z axis,  the  farther  apart  are 
integration  points  along  the  strut.  This  is 
the  reason  for  the  size  restriction  discussed 
under  "limitations." 

• Due  to  the  analytical  modeling  done,  any  part 
of  a strut  whoso  projection  falls  outside  the 
aperture  will  be  ignored. 

• There  must  be  NS  of  these  cards. 
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PHI  CONTROL  CARD 


PURPOSE:  To  allow  the  user  to  specify  which  pattern 

cuts  he  wants. 


CARD: 


/ " 

20 

30 

Ta 

50 

60 

70 

80 

NP 

ANP  (1) 

ANP  (2) 

ANP ( 3 ) 

The  n-.?tr.bprs  along  the  top  refer  to  th?  last  column  In  each  field. 


PARAMETERS:  ( I integer,  F decimal  number) 


NP  (I)  - Absolute  value  is  the  number  of 

phi  cuts.  If  NP  is  positive,  cuts 
will  be  evenly  spaced.  If  NP  is 
negative,  cuts  can  be  arbitrarily 
spaced . 


ANP(l)  (F) 
ANP (2) 

* 

ANP  (7) 


Information  on  phi  cuts:  if  evenly 
spaced  case,  ANP  (1)  is  starting 
phi  value  and  ANP  (2)  is  increment 
between  successive  phi  values.  If 
arbitrary  spaced  case,  the  first 
NP  values  of  ANP  are  the  phi 
values . 


NOTES: 

• No  comments  are  permitted  on  this  one  card 
since  all  possible  fields  are  used. 


• If  arbitrary  spaced  phi's  are  used,  limit  is 

7 different  values.  This  limitation  does  not 
apply  to  evenly  spaced  cases. 
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THETA  CONTROL  CARD 


PURPOSE:  To  set  spacing  and  limits  of  field  values 

within  a phi  cut. 


CARD: 


/ “ 

20 

30 

I 

NTHETA 

THONE 

DTHETA 

The  numbers  along  the  cop  refer  to  the  last  colc.Tn  in  each  field. 


PARTVMETERS:  (I  integer,  F decimal  number) 

NTHETA  (I)  - number  of  theta  values 

THONE  (F)  - first  theta  value 

DTHETA  (F)  - increment  in  theta  values 


NOTES: 

• THONE  should  be  assigned  0 to  obtain  proper 
normalization  of  secondary  pattern. 
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PLOTTING  CONTROL  CARDS 


A special  group  of  cards,  with  a different  mode  of 
operation  are  used  to  control  the  plotting  of  data.  This 
input,  consisting  of  the  AXIS,  TITLE,  CURVE,  and  NEXTPHI 
cards,  can  be  rearranged  to  provide  a wide  variety  of 
plotting  options.  Much  of  the  data  read  in  by  these  cards 
has  an  assumed  default  value.  If  the  user  does  not  provide 
some  of  this  data,  the  program  will  assume  a useful  value. 
Hence,  the  user  has  the  option  of  using  standard  size  plots 
by  entering  only  several  numbers,  or  else,  choosing  his  own 
plot  format  by  entering  additional  data.  Once  the  user  has 
entered  this  data,  it  becomes  the  new  assumed  values  for  the 
rest  of  the  program  and  this  data  does  not  need  to  be  entered 
a second  time  unless  it  is  desired  to  be  changed. 


AXIS  CARD 


PURPOSE: 


To  cause  the  plotter  to  draw  a new  coordinate 
axis  for  the  patterns  to  be  subsequently 
plotted . 


CARD: 


/ 

20 

30 

40 

50 

60 

AXIS  1 

XLENGTH 

YLENGTH 

VLOW 

VHI 

PHILAB 

The  numbers  along  the  top  refer  to  the  last  column  in  each  field. 


PARAMET1:RS:  ( I integer,  F decimal  number) 

AXIS  1 : The  word  AXIS  in  columns  1 thru 

4 and  the  number  1 in  column 
10  signifies  which  type  of  plotting 
card  has  been  read. 


XLETiGTH  (F)  : length  of  the  horizontal  axis  (theta 

angle)  in  inches 


YLENGTH  (F)  : length  of  the  vertical  axis  (field 

pattern)  in  inches 

VLOW  (F)  : value  of  the  field  pattern  in  db 

to  labeled  on  lower  point  of  vertical 
axis . 
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AXIS  CARD  (cont.) 


VHI  (F)  : value  of  field  pattern  in  db  to 

be  labeled  on  highest  point  of 
vertical  axis. 

PHILAB  (F) : decimal  code  that  controls  plotting 
of  the  phi  values.  The  values  are: 

1 = phi  value  for  each  individual 

curve  will  be  plotted.  Use  this 
option  when  plotting  different  phi 
cuts  on  the  same  axis. 

0 = All  the  curves  on  a single  axis 
are  for  the  same  phi  cut  and  phi 
will  be  labeled  only  once  for  the 
graph. 


NOTES : 

• The  default  values  for  this  card  are: 

XLENGTH  = 9.  VLOW  = -72. 

YLENGTH  = 6.  VHI  = 0. 

PHILAB  = 0. 

Thus  the  standard  default  plot  has  an  axis  6"  x 9" 
with  the  6"  vertical  scale  ranging  from  0 down  to 
-72  db.  Only  one  phi  cut  angle  label  is  provided. 

• Labels  on  graphs  make  complete  graph  larger  than  just 
dimension  of  coordinate  system. 


CURVE  CARD 


PURPOSE: 


CARD: 

(curve 


PARAMETERS 

CURVE 

COMP  ( I ) 


To  cause  the  computer  to  plot  a selected 
component  of  the  most  recently  computed 
phi  cut  and  label  it. 


10 


20 


2 COMP 


Ihe  numbers  along  the  top  refer  to  the  last  column  in  each  field. 


( I integer,  F decimal  number) 

2  : The  word  CURVE  in  columns  1 thru 

5 and  the  number  2 in  column  10 
signify  which  type  of  plotting  card 
has  been  read. 

: Component  of  field  that  is  to  be 

plotted.  These  component  numbers  are 
the  same  order  as  the  magnitude  in  the 
output  printout.  The  values  are: 

1 = for  strut  scattering,  principal  component 

2 = for  strut  scattering,  cross-polarized 

component 

3 = for  feed  scattering,  principal  component 

4 = for  feed  scattering,  cross-polarized 

component 

5 = for  reflector  + direct  feed,  principal 

component 
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CURVE  CARD  (cont.) 


6 = for  reflector  + 

component 

7 = for  total  field, 

8 = for  total  field, 


direct  feed,  cross-polarized 

principal  component 
cross-polarized  component 


NOTES: 


• Any  number  of  curves  from  the  same  phi  cut  can  be 
plotted  in  the  same  graph  by  successive  curve  cards. 

• Curves  for  strut  scattering  and  feed  scattering  will  be 
plotted  only  as  far  as  0 = 90  even  if  the  maximum  value 
of  9 is  larger.  This  is  because  these  quantities  are  only 
calculated  in  this  range.  For  6 =»  90",  the  reflector  and 
direct  feed  and  total  fields  are  the  same  curves. 

• Sections  of  curves  that  would  go  beneath  the  minimum 
plotted  values  VLOWspecif ied  on  the  AXIS  card  are 

drawn  as  straight  lines  at  the  bottom  of  the  graph.  This 
situation  occurs  for  very  deep  nulls  and  for  cross-polariz 
curves . 

• AXIS  card  should  precede  the  first  CURVE  card. 


PARAMETERS:  { I integer,  F decimal  number) 


title  4 


NC  (I) 


The  word  TITLE  in  columns  1 thru 
5 and  the  number  4 in  column  10  signify 
that  this  and  the  next  card  provide  title 
information . 

The  number  of  characters  in  the  title 
which  will  be  given  on  2nd  card 


SECOND  card  : 

title;  the  first  NC  characters  in  this  card 
will  appear  centered  beneath  the 
coordinate  axis. 


'•  If  NC  is  specified  as  80  the  entire  second  card  will 
be  the  plotted  title.  The  variable  has  no  default 
value . 

• The  size  and  shape  of  the  reflector  and  the  frequency 
are  automatically  plotted  in  the  upper  right  corner 
by  the  AXIS  card.  The  TITLE  cards  are  useful  to 
record  the  struts  and  feed  details  on  the  plot. 
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NEXTPHI  CARD 


PURPOSE: 


This  card  causes- the  next  phi  cut  to  be  calculated 
after  processing  plot  cards. 


CARD: 


:NEXTPHI 


The  nu.nbers  along  the  top  refer  to  the  last  column  in  each  field. 


L 


PARAMETERS:  ( I integer,  F decimal  number) 

NEXTPHI  3 : The  word  NEXTPHI  in  columns  1 thru 

7 and  the  number  3 in  column  10 
signify  what  type  of  plot  card  this  is 
There  are  no  other  parameters  on  this 
card . 
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NLXTPHI  CARD  (cont.) 


NOTES : 


• This  card  is  necessary  to  specify  that  the  computer 
should  calculate  the  next  phi  cut  it  has  been 
previously  given  and  not  expect  another  plot  command 
to  plot  a curve  for  the  present  phi  cut. 

• There  should  be  as  many  NEXTPHI  cards  as  there  are  phi 
Even  if  no  plotting  is  desired  the  NEXTPHI  card  must 
appear.  The  program  has  been  written  so  that  if  only 
NEXTPHI  plot  cards  are  used  (no  plotting  is  desired) 
then  no  calls  to  plot  routines  will  be  executed  and 
no  output  plot  file  need  be  provided.  The  program  is 
also  written  so  that  if  plotting  is  never  desired, 
only  1 card  need  be  removed  from  the  main  program  to 
eliminate  reading  NEXTPHI  cards. 


After  reading  a NEXTPHI  card,  a new  AXIS  card  needs 
to  be  read  only  if  a new  axis  is  desired.  In  this 

way,  different  phi  cuts  can  be  plotted  on  the  same 
axis . 


cuts . 
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B.  Example  Computer  Run 


The  input  data  for  an  example  computer  run  are  given  in  Table  I. 
This  example  consists  of  three  cases;  the  printed  output  data  for  the 
first  case  is  given  in  Table  II  and  the  plotted  output  is  given  in 
Figures  5,  6 and  7.  This  example  corresponds  to  a 24"  diameter 
reflector  with  F/D  = 1/3.  For  the  first  use  there  are  three  struts 
with  diameter  2ap=0.37",  3=68°  and  <fi5=90°,  210°,  330°,  respectively. 

A circular  flat  plate  with  radius  c=1.2"  is  used  to  model  the  feed. 

The  pattern  plane  is  <()=-15°  and  the  frequency  is  11  GHz. 

The  printout  of  the  first  case,  given  in  Table  II,  first  echoes 
the  input  data,  and  then  prints  an  interpolated  input  pattern  that  has 
been  deduced  from  the  arbitrarily  spaced  primary  feed  pattern  data 
points.  The  printout  of  the  coefficients  of  the  scattered  wave  was 
used  in  the  development  of  the  program  and  can  usually  be  ignored  by 
the  user.  The  bulk  of  the  printout  is  the  different  contributions  to 
the  secondary  field  of  the  reflector.  In  the  sample  a theta  spacing 
of  1/4°  was  used  so  that  the  corresponding  plot  would  have  very 
smooth  curves;  a wider  spacing  can  be  used  with  increased  roughness 
of  plotted  curves.  The  principal  and  cross-polarized  components  of 
each  contribution  are  referenced  to  the  Huygen  polarization  reference. 
After  90°,  only  the  reflector  and  direct  feed,  and  total  fields  are 
computed  and  listed.  The  printout  for  the  rest  of  the  three  cases  is 
similar  and  has  been  omitted. 

Three  examples  of  the  computer  plotting  subroutines  are  shown  in 
Figures  5,  6 and  7.  The  example  in  Figure  5 demonstrates  the  capa- 
bility to  plot  the  individual  scattered  field  components:  the  basic 

contribution  from  the  reflector  and  direct  feed  radiation,  the  strut 
scattering,  the  feed  horn  model  scattering,  and  the  total  field 
which  includes  the  various  scattered  field  components.  In  this 
example  only  the  principal  polarization  components  are  shown.  The 
antenna  geometry,  pattern  plane,  and  frequency  correspond  to  the 
example  computer  run  listed  above.  In  this  plot,  the  standard  size 
and  scale  have  been  used. 

The  NEXTPHI  3 card  causes  more  computations  to  be  performed. 

Since  only  one  value  of  phi  has  been  specified  above  and  since  NCASE 
was  read  as  3 the  next  card  is  read  as  specifying  a new  strut  and 
blockage  problem  with  the  same  reflector  geometry  and  primary  feed. 

A new  strut  configuration  with  a sinqle  strut  at  phi=90°  has  been 
set.  The  second  case.  Figure  6 shows  the  capability  to  plot  dif- 
ferent pattern  plane  cuts:  (i)=-15°,  0°,  90°  are  shown.  In  this  example 

only  the  principally  polarized  component  of  the  scattered  field  from 
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the  struts  are  shown.  Again  the  antenna  geometry  and  frequency 
corresponds  to  the  above  computer  run.  Each  NEXTPHI  3 card  causes 
the  phi  values  defined  on  the  phi  value  card  to  be  executed. 

The  capability  to  calculate  cross-polarized  components  is  shown 
in  Figure  7 where  only  the  total  field  is  shown  in  both  principal 
and  cross  polarizations.  The  reflector  geometry,  the  direct  feed 
pattern  and  the  frequency  are  the  same  as  above.  However,  instead  of 
tripod  struts,  there  is  only  one  strut  located  at  ip^=20°.  The  strut 
diameter  is  2.0  inches  (1 .86X ) and  6=68°.  The  strut  length  is  the 
same  as  above.  The  strut  angle  4>s=20°  and  the  large  strut  diameter 
was  chosen  to  demonstrate  a large  cross  polarized  component  of  the 
strut  scatter.  This  cross  polarization  effect  dominates  for  <)!>45° 
in  the  (}i=120°  plane  which  is  near  the  scattering  cone  of  the  strut. 
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TABLE  I - Input  Data  Card  Listing 
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TABLE  I - Continued 
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TABLE  I - Continued 
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TABLE  II 

COMPUTER  OUTPUT  LISTING 
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C.  Appendix  - Computer  Program  Listing 

The  listing  for  the  computer  program  described  in  this  report  is 
ven  below: 

"Circular  Reflector  Antenna  with  Feed  and  Strut  Scatter"  (OSUPATT) 
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♦ NUMPfR  {'if  LlNfAR  PFED  CPfFFlCIFNTS  OSU 

O'.U 
f «;  j 

os  SfiFI  niFFRACTION  CnFFFlCIfNT  CSlI 

OM  JURO  oiffraction  coefftcunt  nsu 

osu 

KLtWA.IlOF  : RfF»R  TO  SUBROUTINF  GTO  HSU 

f'StJ 

f’Ml 

DIMENSION  (.(  lOOJ  ,hMOO)  , ZY  MOO)  » ;x  f 100  ) ,DLF  < SOI  ,OLH(  SO)  ,HF  T ( lOD)  , hSn 
i)  F MOO  1 ,0M  1 (>0  I.  TP(  1 00  I . TM(  100  ) tFXF  ( SO  I ,PXH(  SO  I , FFF  (SOI  ,FFH(  Sni  . OSO 
t B J(  ?S  I , ANP  (l0),tHM(?l,MS?(A),ZS?(^.),RLn(4|,PHlS(4|,PLnTST(fl,7?l|  OSU 
CDMPL  f X J ,CP  ,CP1  fCP?»i)S  I »nS?,OHl  .OH?,F  01  ,FP7,FM  ,FT?  *FFT  ,FFP  ,FF  , OS II 
1C  P*.  ,0S  , ,0HC  ,A  1 ,A?,A^,Bl*B?,0‘^,FY,Ff,FP,CP‘^tAM|?S)  , HNK  ( ?S  | , OS  U 

lHKP,F>^,f  Z,l  SY,fST,fSP,ffHF,FPfU*FrT,FTP,CPP,nM(?Sl,FAT,FAP*CPOt  OSU 

IF  f BY  ,M  p T ,l  M*P  .1  SM  ,(  STP  ,f  PL,MPL  HSU 

KF  Al  MAX  DSl) 

COMMON  /!  TM  / J.CPO,  AM  ,PM,P  I ,RK  ,MX  ,DFL  f A , TF,  OF  ,Rf  TN  (\SU 

FORMA  r ( IF  1 0 .4  , [ M nsu 

fOKMAF  ( //•  riKCUlAR  RFFLFCTOR  ANTENNA  FAR  FIELD  PATTERN  WITH  OSO 

fSCAlTFHiNG  FROM  f(FD  STRUTS  AND  F F F 0 * , // 1 S ♦ * AN  T F NN  A OF  OMF  1 R V * , //  OSlJ 
IT  10,  • »F'F  RTUMI  OIAMITFR  =*  , F H . , TS  0 , • L AMOA  , F .1 . 3 , T 7 S , • 1 NC  H F S ' , /DV  1 1 

1T?4»,F  M.  T,  y.y,F  P.  T zs  , *MF  TF  RS*  r/T?«V,Fn,3  , T 7S,  •WAVELENGTH*  ,// no,  r»S  U 
1'FOCal  distance  to  niAMriER  ratio  ^ , FR.  1, //T  lO.  • fpe  QUENCY  s»,  osii 
IFK.A,'  MHZ  • , //T  10, 'NUMfitR  OF  APFRTURF  POINTS  s*,!S,///|  OS  U 

4 FORMAT  n?M  *STRIJ1  SC  A T • , T SO  , • F F E D SCAT»,T7?,»  rffi*dirfct  fffpsu 

HU  •,Tno,*iniAi  Fnin*/  vx  ,a(  ox , * pr  tnc  • ,n  x , »cpnss  • ,?x  1/  psu 

X ox , n ( 4 X , 'DU*  ,sx , *01 c M I nsu 

OSlI 

READ  ( *^ , 1 ,F  N(>3  10  1 I nZ.F  REO,FDR  ,NAP  OSIJ 

tlSU 

IF  (NAP. n. ‘01  NAPsSO  OSl» 

IF  (NAP.CT  NAPr^O  flMl 

p|  M.  nsn 

o^il 

j- (0. , 1 , 1 nsu 

M A X -U  . OSlI 

Ri AMs  ADO . /F  R r 0 OSU 

FMi('.(*?Ss  nsu 

K(  AMl - «L  AM/r  M OSU 

01 -0?*FM  OSU 

psHl/Rl AM  OSU 

WR  n F (S  , U?  ,R  L A«1  ,n  I ,RL  AM  ,f),  FDR  ,FRE  0 ,NAP  OS'* 

A-!)/.*.  OSU 

E r f OR  #fi  OSU 

Ml  rNAP«  1 flSU 

0!  I i A /NA  P U 

P = r«A*A/(A.*M  OSU 

ROiSuK  I { A »A  ♦ P«*B1  0<  H 

( pu=c I xp ( -J *RK *R ni  USU 

Cp=J*CP(I  OSU 

CPssp|*A*CFxp(j*(pi/4,-HK*Rrn)  nsu 

Al  ps  A T AN.  ( A , P I F»SU 

alpha  Al  P*lf.C  ./P  I nsu 

V INA  s IN ( Al  f I ' Su 

Co*.  A-C(l‘  ( AL  p I (■'  u 


• IINEAR  FEED  CntFFlCHNTS  ♦ 


OSU 

10 

nsu 

?0 

f'SU 

30 

OSU 

40 

OSU 

SO 

f S J 

CO 

csu 

10 

nsu 

PO 

OSU 

VO 

nsu 

100 

OSU 

1 10 

f’SU 

1 ?0 

usu 

1 30 

0S(/ 

t 40 

OSU 

l‘,n 

nsu 

1 ^0 

os  u 

1 70 

nsu 

1 HO 

nsu 

100 

OSU 

?oo 

(\S  u 

710 

nsu 

?20 

nsu 

?30 

nsu 

2rO 

/nsu 

?so 

ns  u 

?co 

usil 

? FO 

OS  u 

?flO 

FOSU 

?vo 

OSU 

300 

OSU 

3 to 

OSU 

3?0 

OSU 

330 

(ISU 

340 

OSl» 

3 SO 

fi‘  u 

3zO 

nsu 

370 

nsu 

3H0 

nsu 

■»vn 

nsu 

4 00 

OSU 

410 

nsu 

4?0 

OSU 

<♦  1O 

OSU 

440 

nsu 

<.*■0 

OSU 

4 1 n 

OSU 

4 fO 

nsu 

A«0 

fjSu 

400 

M 

SOO 

OSU 

S|0 

n<  u 

*c 

(iSU 

S30 

OSU 

».40 

nsu 

4S0 

nsu 

s^O 

OSU 

S70 

'SU 

M«D 

(«'  u 

* » 0 

mu 

Z/00 

f.SU 

t IC 

rsu 

/ .’0 

OSU 

t.  50 

• • ■ , t,- 

wLl: 


rior*oo 


C 

C 

10 

1 1 


c 

c 

c 

c 

c 

1 


16 
1 7 


C 

C 

C 

C 

c 


70 

71 


7? 


P F 40  I •»  t 1 r 1 N7 
f OkMA  T ( I 10) 
wP  I n < 6,  U ) N7 

format  iMOt'NUMprn  nr  input  fffo  valufs  -**I7,//) 
call  LNMHOl  F *P>0  .Ff  f ,N7) 

C4fL  LNFnn»LM,PXM,FFH,N7) 

N3  ^N7-l 


* FFFO  PATTFPN  ♦ 


WRITF  16,  n ) 

F0»^MAT  (//110,*FFFn  PATTFFNS//T?^,*PS1*,T37,*FF%T5  1,»FF0B*»T66, 
$ •F^^•  , TOO  , M Hf>P  * , //) 

PS1=0- 

On  IF  L=  1 ♦ ’7 

call  L sum (PS  1 ,PXF,FFF,0LF*N3,FF) 

CALL  LSUM(PS1  ,PXH,F  FH,0LH,N3,FH) 

IF  (FF.FC.O.)  FM)D^-600. 

IF  (FH.FQ.O.)  F hop =-500. 

IF  (FF.Nfc.O.)  FF0n=?0,*ALnG10(APS(FF)  ) 

IF  (FH.NF.O.)  FMOE=70.*AL0010( ABS(FH) ) 

WRITF  16, \6)  PSl  f FF  ,FF0R  ,FV(  ,FV^^P 
FORMA  T ( 1 IP , FR.?  ,AF 15.A  ) 

PSI=PSI*5. 

CONT I NUF 

• niFFRACTION  J 

GAM  = A T AN? { 7 . AF  ,A  ) 

GA««A  =GAM*1 PO. /P  I 

PmIPs  1 P0.-(-A“MA-ALPHA 

5=100. ♦Rn**7 

SL=RO 

1L0P‘P 

w A=0.0 

PvnP0='^60.-0A”^*A 

CALL  GTLM  PHI  PO,PMIP,QO.,^,PO,WA  , 5L  ,AS  , HOP  ,nS,OH) 

AN  Is- 1 PC.  - AL  P»^A 
an7- I PO .-gamma 

• fffo  sc ATTFRING  MOOFL  • 

• 4 444^*4***««4««4«4**«**« 

PFAf)  (5,71)  NCASl 
IF  ( N(  AS  * .LI.  1 ) NCASr  S 1 
on  I 00  NC - 1 ,NC A?  F 
WRlTf  (6,7(M 

FORMAT  l//T10,*Frf()  SCATTFRING  MOOFL*,//) 

RfAn(5,71)  KG,  in  IM ( I ) , 1 = 1 ,KG  ) 

F ('RMA  T 1 I |0 , 7F  1 0 .5  ) 

o(  )=urM  n -M' 

Of  =0C  I /RL  am 
C - DC  / 7 . 

IF  |K(, . F C . I ) GO  TO  7 T 
AV-rnC 

P.  Y I r f)  J M ( / ) ♦ F M 
P V - p V I /F’  L AM 

f PPr  A X ♦ H Y *(,0  *r  F X P C - J ♦R*'  ♦!»  ) 

WPITI  16,77)  DIM  n > ,T')M1  7 ) ,DC1  , PYl  , AX  ,PV 

F nu*1A  T ( 1 I ' , • W F I I AN(.UL  AR  FIFO  HLOTKAf.l  »,//T^0,  ‘WIDTH  r‘,FP.'»,T*‘S, 
1'hMGMT  : • , f p.  >,  T60,  ' I Mf  H»  s * , / T 77  , *■  n . T , f ^ ^ . r M . 1 , I no  , • MFTIP**,/ 


nsn 

6R0 

osn 

650 

nsu 

660 

mn 

670 

PSU 

f PO 

f'CU 

6 90 

nsn 

700 

nsu 

710 

nsn 

770 

nsu 

730 

nsu 

7 AO 

nsu 

760 

nsu 

7 60 

nsu 

770 

nsu 

7P0 

nsu 

790 

nsu 

POO 

nsu 

PIO 

nsu 

P70 

nsu 

p'»0 

^su 

p/.O 

usu 

P60 

05U 

P60 

nsu 

P70 

nsu 

PRO 

nsu 

R90 

nsu 

900 

nsu 

9 10 

nsu 

970 

nsu 

930 

nsu 

94.0 

nsu 

960 

nsu 

960 

nsu 

970 

nsu 

9P0 

nsu 

990 

nsu 

1000 

nsu 

1010 

f^<U 

1070 

nsu 

1070 

nsu 

1040 

nsu 

10*^0 

nsu 

1060 

nsu 

1070 

nsu 

ICPO 

nsu 

10  90 

nsu 

1 1 on 

nsu 

1110 

nsu 

1 1 70 

nsu 

1 1 30 

nsu 

1 140 

nsu 

1 1 60 

nsu 

1 1 60 

nsu 

I 1 70 

nsu 

1 ] «0 

nsu 

1190 

nsu 

1700 

nsu 

1710 

nsu 

1770 

f SU 

1730 

nsu 

1 740 

nsu 

IP^-O 

f^SU 

1760 

OSU 

17  70 

nsu 

1 ?PO 

\ 1 

1 790 

n‘  u 

1 ^00 

tiv  u 

1 MO 

80 


c 

r 

c. 

c 

c 

31 


IS 

C 

c 

c 

C 

c 


A1 


S’ 


tl  w,(  f\.  3 , Ih  < ,f  H.  IHO,  ' w*Vf  LI-N(;3H*  ,// I nsu 

&n  TO  ?s  n'.i) 

wRlif  (S.7SI  n t«  { I I ,f)C  1 .nc 

MIRMAI  ( 1 1 S , •(  IR  CUl  » R flfn  BtOCK  »&r  > , //T?S, '01  »MFTFR  -',rn.3,T5S,  o^u 
»•  INUir  5,  • ,/1  3F  ,c  ?.  5 , TSS,  • Ml  tFRS*,/T3B,FB.3,FSS,  • WA  VFL  F N'VTM  • , // 1 PMI 
Ci’l'r7.*PI»r*»7«CI'(IM-J»RK*H)‘CP  ns  II 

coNTiNur  nsu 

nsu 

4. •««**««»*********•«***  nsu 

• I tin  STRUT  Lnr.ATioNS  * nsu 

444*44««444*4*44*4**4**4  HSU 

nsu 

RFAf)  (s,31)  f)P;,PFTA,NS  nSII 

FORMAT  I'HO.S.IS)  nsu 

nPirOP7*FM  nsu 

np-ftp  1 /Ri  4»t  nsu 

AP-OP/7.  nsu 

WRITF  (6,37)  nP?  ,RF  T A ,NS  , DPT  ,0P  OSU 

Fnp«4I  (110, •FIFO  STRUT S • , 7/T7 A , *0 I4MFTFR  -'.Fn.T,'  INCHFS',T6n,  OSU 
t'FiFTA  - •,i(..i,'  nrr,RF  F s • , TOO,  • numbfr  of  struts  =•  • ,i?,/t3a,fb.3,  nsu 

»•  MF  IFRS  ' ,/T  3S  ,F  fl.  A,  • WAVFirNl.TM  • ,// 1 nSU 

on  3s  T=i,NS  nsu 

RfAn  (S,3A)  PS77,7S77,PL07,PHIS7  OSU 

FORMAT  (AFl(l.S)  nsu 

RS?1=RS77*I M nsu 

7S71’7S2?*FM  nsu 

RLin  -RLO?*l  M nsu 

RS7( I )=PS71/Rl  AM  nsu 

7S?(I  )=^7S7I/RLAM  nsu 

RLO(  1 I =R  LOl  /RL  AM  nsu 

WP.ITF  (6,3A)  I .PUIS7  ,RS?7, 7S77  ,RLn7,RS?l  . 7S71  ,RLni  ,RS7I  I 1 ,7S2(  T 1 , nsU 
♦ P(  n(  I ) nsu 

FORMAT  (T7S, 'STRUT  NUMB  f R • , 1 7 , / / T S 0 , • PHI  S =',F5.I,'  PFC-R  FF  S ' , T B 0 , OSU 
»•  RHUS  ',T<J3,'7S',T103,'LFNr.  rM',//T73,3F17.3,TllS,'INCHES',/T73,  OSU 
»3FT2.3,T11S,'MFTIRS',/T73,3F17.3,T1 1 S , • WA VEL FNGTh • , // I OSU 

PMisi  I ) = PHi s7*Pi /I BO . nsu 

continue  nsu 

nsu 

• coefficients  of  cylinder  scattfrino  * nsu 

nsu 

PF  tn’BETa«pi/i  no  . nsu 

SINB^SINIBE IN)  nsu 

arc-Rk»ap4SInb  nsu 

WRITE  (6, A! I ARG  nsu 

FORMAT  ( //T 1 O, 'A PG- • , E B .A , //)  OSU 

Mx=^ARG»  1 1 nsu 

nn  A3  M = i ,Mx  f'Su 

L’M-i  nsu 

CAIL  Bf s J ( ARC , I , BS J ,0.001 , I ER ) OSU 

call  ru  sy  ( apg, l , bn , I f r ) nsu 

BJ(M|-RSJ  (>SU 

)'NK(M)  = ns  j-j-*BN  nsu 

AM  ( M ) =-H  J ( M ) /ONK  ( M ) OSU 

I E (M  .0  ) . 1 ) GC  TO  A7  nsu 

AM  1 1 ) -O.S*AM(  1 ) CSU 

f,n  TO  A3  OSU 

B JP’H  J(  L )-(  L/ARG  )»I1J(  M)  nsu 

MI'P=IINY  ( L )-  I I /arc.)  ‘OMR  ( M ) nsu 

PM  ( M I I-B  JP/HF  P nsu 

IF  (arg.gt.1.)  gfi  in  A’  nsu 

P PA=R|  A(  ( AM  ( M)  )/RF  Al  ( AMI  1 I ) FSU 

IF  (PR  A .(.T  . 1 . OE-B  ) GO  TO  A3  OS  U 

PI  A = A I MA  (.(  AM  ( M 1 ) /A  1"AG(  AM(  I ) ) OSU 

IF  (P  I A .f  T . 1 .OF-F)  F.n  TO  A3  OSU 

GO  H'  AA  T'SU 

CONTI NUf  OMI 


1370 
1 3’0 

I3A0 

I 360 
1 3 60 
13  70 
1 3B0 
1 TOO 
lAOO 
I A 10 
1 A70 

IA30 
lAAO 
lASO 
1A60 
IA70 
1 APO 
lAPO 

I son 
isin 
1 S70 

15  30 
ISAn 
1S‘  0 
1S60 
1670 
ISBO 

1 6 on 
1600 
1610 
1670 
16  30 
16A0 
I6sn 
1660 
16  ro 

1 6M0 
16P0 
I 700 
1710 
1 720 
1730 
I 7A0 
I7'0 
1 7 60 
1770 
1 7B0 
I 700 
1 BOO 
1 filO 
I B20 
IB  30 
IRAO 

1 nso 

I B60 
I A 70 
1 RBO 
Ipoo 
loco 
I'MO 
1020 
10  30 
lOAO 

I'76f) 
1060 
10  70 
|OB0 

1 oon 


r r- 

•iJlt 


copy 
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“>r>onr> 


C 

C 


ss 


^ I 


C 

C 

C 

C 

C 


BM \ \ \ = -0 .^♦R  ,M  ?) /MNK  t ? ) 

w«?  m 

hOH‘^/V  T ( /I  h,  *(■  Of  F F Tf.  I FNTX  OF  IHF  <;rftTTFRFn  W4VF  FROM  A CYLINOFF', 
t//l?0, *M', ’AMIM)* ,T6S, 'RMfM) *,//) 

WK  nr  1 S ,‘.F> ) U1  M > ,BM  1M  > = l ,MX  ) 

f nWMAT  <1  W. , 1 ,4  f I 5 ) 


RFAn(*>,‘3)  NP,  UNPfK  >*K  = 1 ,7  > 

f OKMAT  ( ! lOt  ?F  10  .**) 

«fAn(‘**sn  fnHfTA,  thone  * othf  ta 
I f vr^=  1 

w .Ll.  0 1 ifvfm=o 

NPr  I A { NP  1 
on  K=  1 ,NP 
PMIsANP  (♦•  i 

If  mvfu  .10.  n poi  = ANPm»iK-n*ANPi?) 

WR|T»  (6,Sf)  PHI 

fnw^^AT  ( // ,no**  PHi=  • ,FF..? , • nFGPFfS*,//) 
PmT  =Ph  I •»  P I / 1 HO  . 
s 1NP=S  U4  1 PMU 
corp-cns(PHi ) 

C0S?P  = COS (? .♦PHI  ) 
p-S1N(?.^phI  ) 
s.  1 N4P=  ^ I N W*  . ♦ pvn  1 
cnr,/.prcnr.<A.»pHi  > 
FAL=^SO«ni.-ISlNA*MNP)A*?l 
F 1 A=COSA  **j  I NP/»  A L 
FPA^CnSP/FAL 


♦ APFRTORF  FIFLOS  • 


RMO-0 . 

on  M I M , M I 

ppTf  -RMO*RMO/f  ♦F  ) 

R P = SOP  1 ( RHn^RHfU  PP ♦RB ) 

$ 1 =A 1 AN? I RMO,BB) 

PSI-S1*1P0./‘’1 
C0SS-=CCS1  \ 

* INS  = S1N(ST  ) 

R»  Til  )r<;t.’RTI<l.-(<;iNS^SlNP)<*?)) 

Tp(i j-i.«crsr 

TM(  I ) = \ .-C(»SF 

CALI.  LSO‘-|PSI  ,P»fSFrM,niH,Ol,FHl 
CALL  LMJM  {PS  I , PXF  ,FF  F , DLF  *N3,FP  » 
CALL  LSU^*  f P'  I t PXH,  r f H,  f)L  3,FH) 

OF  ( n f Ff  -F>^  ) *F/  ( ?.  *RP  ) 
PHOrRHO»(>F  L 
CONT I NUf 
TS=TF (MI ) 

01 =0F ( YI  ) 

PP-TS-OL  *CnS?p 
WRITF 


♦ PATTfPN  CALCOLATJON  • 


no  on  N- 1 ,\THF  TA 

THFTA  = TMO{j1  ♦(N‘1  1*(^MFTA 
TMF - TMf  T A*P I /I  PO  . 
MNl-StNi  I 
cnsT^cns ( f Hf  1 

( p i^crxpf  .i*RK4p4fnx7  ) 


rsii  ?roo 
nsu  rn 10 
nso  ?o?o 

fiSO  70^0 

nso  ?r/«n 

fP.O  ?0‘.0 

nsu  ?c^,o 
nsu  7070 
OSU  70P0 
fiSU  7(»90 
nsu  7100 
nsu  7110 

nsu  7170 

f'Sif  71  xn 

nsu  71  AO 
nsu  7)*»o 
nsu  7160 
nsu  7170 

nsu  71H0 
osii  7 1 on 
nsu  7700 
nsu  7710 
nsu  7770 
nsu  7730 
nsu  77AO 
nsu  7750 
osu  7760 
nsu  7770 
nsu  77PO 
nsu  ?7«*n 
nsu  7300 
nsu  7310 
nsu  7370 
nsu  7330 
nsu  73AO 
nsu  73'^0 
nsu  7360 
nsu  7370 
nsu  73RO 
nsu  7‘»90 
nsu  7A0O 
nsu  7A10 
osu  ?A?0 
nsu  ?A?o 
nsu  7AA0 
nsu  7AS0 
nsu  7A60 
nsu  ?A70 
nsu  7AP0 
nsu  7A90 
nsu  7500 
f!SU  7610 
osu  7670 
nsu  rs‘*r 
nsvj  7640 
nsu  7650 
fisu  7 6 60 
r»U  7S70 
PSU  7*^«0 
nsu  7690 

nsu  7600 
nsu  7610 
r*^u  ?6?o 
rnu  7630 
nsu  76A0 
u 7650 
n^u  7660 
osu  76  FO 


82 


loonooo 


t 

C 

C 

C 

C. 


6' 


66 


70 


71 

77 


r 


n nHF7».iT.op.  I GO  in  6? 

IF  (X.LT.?.*>t  r.r  TO  76 
GO  TO  70 

IF  I » .GT  . 74  . ) GO  TO  70 


» 6[’irturt  I NT  fgr  »t  inN^» 

RHO=0. 

no  6G  ]=  1 .“FI 

xr-RK*RHO»FINT 

If  (XP.GT  .4.rT-7  I GO  TO  6X 

V ( XP*XP  ) 74 . 

M JO'l .-Y«Y*Y/4 . 

BJ7=Y/7. 

0,14x0. 

GO  TO  64 

C Al  L Of  S J (XP  ,0  .P.KsO.OOl  , I FR  1 ) 

CALL  Of S,MXP,7,0 J7.0.00I . IFRl) 
f,  AU  0ESJ(XP,4,f\J4,0.nOl  .IFRII 

0(  I ) = -($! N7P*TI  ( I l♦PJ7♦n.G»SIN4P♦nFI n*Rj4)*RH0*TM(l  )/t7.*BET(] ) I 
G(ii--nFni*pjo*OF(i)*nj?»cn5  7P)*RHO 
RtnixPHOxPFL 
CONTI  NUT 

CALL  OSF  ( OF  L ,H  ,7  X.M]  ) 

Call  OSF  (Of  L ,G,7Y,pI  I 
F XxCP*RK»COST • ZX (M I I 
t Y-CP*RX*C0ST»  ZY  (MI  ) 

FZ--CP«RX»SINT*(C0S'>*7X(MII«SINP*ZY(M1)T 
f n = _f  X*  S I NP  »f  Y»COSP 

rixFX»COST»cnSP»FY*COST»SINP-F2»SINT 
U I TMFT A .Nt .0 .0  I GO  TO  77 
MAXr70.*ALOGIO(CAF(S(FYn 
WRITF  (6,66)  max 

format  ( //t 16, 'AXIAL  fiflo  =',fr.?,'  no‘,//) 

GO  TO  77 


* TWO  POINT  MFTHOn  * 


T7  = S0R  T ( A/SINT ) 

PS  1 X 1 HP, -TM'  1 A 
S1-P1-TH( 

SINS=SIN(SI  I 
COSS=COS (SI ) 

FAr,xSCRT(I.-(SINS»SINP|»*7) 

IF  ( F AC . FO.O . ) CO  70  71 
F TrCOSS»S  INP/f  AC 
FPxcnsp/f  At 
CO  TO  77 
F I = I . 

FP=0. 

CF’i=crxp(.i*Ri<»(A«siNT-Rn)) 

CP7xCFXP(-J*HKa(A*SINT'R0II 

PHI  1= 1 no. -gamma  » THI  I A 

PH  I 7= 1 PO . -GA“MA- TH(  ta 

IF  (PHI  7 ,LF  .n.l.)  PHI  7=  T60  , 'PHI  7 

call  gto(phi  I ,PHiP,tin.,s.Ro,wA,SL,Asi , ilop  ,nsi  ,nHi  I 

call  GTO ( PHI  7 , PH  IP ,P0. , S ,R0, WA , Sl, AS  I , I LOP ,0S7 ,nH? ) 

FTIx-T7*CPI*0MI»PP*FTA 

FT  ? = -.)»T7»CP7*nH?«HP»F  lA 

F p I xt;  »cp  1 ATiSi  *H  p*F  r A 

fP7=^J*T?*LP7*OF?»PP*FPA 


OS  (I  76  HO 
nso  7660 
050  7 700 
OSO  7710 
OSH  7770 
OSO  77TO 
riSO  7F40 
OSO  77SO 
OSO  7760 
OSO  7770 
OSO  77H0 
OSO  7760 

nso  7R00 
osij  7HI0 
nso  7P70 
050  7PT0 
OSO  7H40 

nso  ?Bsp 

OSO  7H60 

nso  7H70 

OSO  7HR0 

nso  7P90 

OSo  7900 
0S(|  7910 

nso  7970 
OSO  79?0 

nso  7940 
nso  7960 
nso  7960 
OSO  7970 
nso  79H0 
nso  7990 
nso  7000 
nso  7010 

OSO  3070 

nso  3030 
OSO  3040 

nso  3060 

OSO  3060 
OSO  30  TO 
OSO  30H0 
OSO  3090 
PSO  3100 

nso  7110 
nso  31  70 

nso  3130 
OSO  3140 

nso  7 160 

(ISO  3160 

nso  7170 
OSO  31  HO 
OSO  3190 
OSO  7700 
OSO  3710 
OSO  3770 

nso  37X0 
nso  3740 

nso  3760 
OSO  3760 
nso  3770 
nso  37HO 

OSO  3790 
OSO  7xno 
OSO  3310 
nso  3370 
nso  33  30 
PSO  7X40 
O' LI  XX  60 


C 

C 

c 

c 

c 


7t. 


7‘ 


C 

c 

r 

c 

r. 

n 


7 7 


• niRrcT  FfFf)  • 


f F=^0. 

}>■  I ThM  A .r.T  . ANI  I f.n  T(l  74 
CALI  LSUM lPSI,PXF,IFF,nLF,N3*Fr) 

C A!  L L SUM  (PS  I ,PX  H,  f FM,  rn  U ,N3  ,FM  ) 
r F • ( ( Ff  ♦ F H)/;*.-CnS?P^  ( F L-FH  )/?  . ) ♦F  *CP3 
t M = -F  F*M 

F F P r F F ♦ F P 

IF  (IHMA.LF  .‘^O.  .n«.TH»  TA.GT.AN?  ) GO  TO  7S 
f P?  = 0. 
r T?=0. 

FP  = F PI  4f  F^->4F  FP 
i T=f  T i + rr  2*t  F T 
G(»  rn  rr 


* RING  CURRfNT  MFTHOD  ♦ 


cont  I Nirr 

CALL  Bf S J( X ,0, BJO.O.OO? , TFR ) 
CALL  0FSJ (X ,2,BJ?,0.001 , rCR ) 
CALL  RFSJ(X,4,BJ4,0.00I ,TfR> 
osc-(7s*cnsT 

A1=TSA (DH-OSC )4DI ♦ (OH ♦DSC )/?. 

A?  -.-TS*  ( OH  4 0 SC  )-0L  ♦lOM-DSC  ) 

A3=nL*(OH40SC)/?. 

OHC=DH*CnST 

Bi  ^TS*(os-OHC)-OLAtr)S4nHC)/?. 
B;*--  f S'-  ( 0S4  OhC  ) 4 OL  ♦(  DJ  -Due  ) 
e3--0LA(0S4nMC )/?. 

fPL-“CP4^(Al#pj04A?»RJ?4A3^HJ4> 
MPL  = CP4*(HIAHjn*B?*F^j?4bi^BJ^) 
|V--rPL»5.  IN  P**?4  HP  L*  CO  SP^^? 
t T— FYASINP 
F P=FY*COSP 
CnNT INUF 


• FMD  F^inCKAGF  SCATTFRING  * 


IF  ( IMF T A .L F .90.  I Cn  TO  RO 
I F PY-0  . 

F ri =F  T 
e T p=f  p 
GO  in  HS 
BO  IF  (KG.FU.l ) on  in  R1 
XA=kK*AX*SINT*CnSP/2. 
Xh=RK*RY^SlNT*SlNP/?. 

IF  (XA.FU.0.0)  GA=1. 

IF  (XA.G 1.0.0)  GA=S1N(XA)/XA 
IF  ( XH  .FO  .0  .(■ ) r,  Br  1 . 

IF  (XB.Cl.O.OI  GB  = SI0(XB)/XPI 
F FPY  = CPP  •(  P iACAAGF^ 

Cu  TO  H? 

PI  r.ONTlNlK 

XC=PK4tA*  UjT 

CALL  ((F  SJ  ( XC  , ! ,e  Jl  ,0.001 , 1 F R 1 ) 
ir  (XC.fO.O.)  F F BY=CPP*CP3/?. 

IF  ( XC  .r.T  .0  . ) f F hY  = CPP*CP3*RJl  / xc 
ft?  Contihuf 

FFf.p^F  FHY^CnSP 
FF  BTsCFRY**  IfiP 


OS  Vi 
OSU 
OSO 
nSu 
nso 

OUf 

nsu 

nsii 

nsu 

nsu 

f>S  0 

osu 

nsu 

nsii 

nso 

osu 

nsu 

nsu 

nsii 

DSU 

nsu 
nsu 
osu 
nsu 
nsu 
nsu 
nsu 
nsu 
nsu 
ns  u 
nsir 
nsu 
nsu 
nsu 
nsu 
osu 
nsu 
nsu 
nsu 
nsii 
osu 
nsu 
osu 
osu 
nsu 
nsu 
nsu 
nsu 
OSli 
nsu 
nsu 
osu 
nsu 
nsu 
ns  LI 
nsu 
nsu 
nsu 
nsu 
iisu 
nsv» 
nsu 
osu 
evsu 

l''SU 

nsu 

nsu 

nsu 


33*. 0 
3370 
33RO 
3390 
■<4  00 
34  10 
34  ?0 
34  30 
34  40 
34SO 
3 4 60 
34  70 
34H0 
3490 
3*^00 
3S10 
3S?0 
3S30 
3S40 
3S‘^0 
3 6 0 
3S70 
3SfO 
3'  on 
3600 
3610 
*>6  70 
?6'40 
3F  40 

3660 
366  0 

36  70 
36F40 

37  90 
3700 
37)0 
3770 
3730 
3740 
3760 
3760 
37  70 
37H0 
3 7 90 
3H00 
3fl  10 
3P70 
3R  30 
3040 
3P  60 
3M60 
3R70 
3FiH0 
3B90 
3900 
3 9)0 
3070 
39>0 
3940 
3940 
3960 
39  70 
39pO 
39^0 
4000 
4010 
40?0 
4030 


84 


(' 

f 

(. 

c 

r. 


81 

80 

OC 


o', 

ICiO 

lOI 


c 


c 


♦ ST«(yT 


»»•******«** 

'CATTf  K7NG  * 


rs  TT= ( 0.  , 0.  ) 

F S TP=  ( 0 . ,0. ) 

on  (10  I I ,rjS 

PtiSrPH  I S ( 1 ) 

7N7=?S?1 1 ) 

RN?''K  S7  ( I I 
PNLI=8in(  I ) 

r tl  L ST8U  I t PHI  ,THFIA,P((S,7N7,RN?.PNn,eST,ESP| 
FF,  IT  = F5TT*(  ST 
FSTP=FST  P«F  SP 
CUNT  irniF 

F n-F  T ♦( ' n *FF(iT 
f TP=FP*F  SIP«(  Fpp 

C Al  L mi  Y SRC  ( FT  PT  ,f  FRP.COf  P.STNP  1 
nOPHS  ( ( f 'IP,  MAX  I 
npPHS(FFBI ,MAX) 
miYSPC  (I  STT  ,FSTP,CnSP,SINP| 

DPPHS ( F STP, MAX ) 

OBPrlS(f  'll, MAX) 

PI  r'TSI  ( 1 ,NI=RF  At  (FSTT) 

PLOIST (? ,N) =Rf  AL  (F  S IP) 

PLCiIST  (3  ,N)  =RF  AL  (FFPT  ) 

PI  nisi  (0  ,N)=RF  AL  IF  FPP  ) 

COMTImif 

CALL  HUYORC  tFT,FP,CflSP,SINP) 

CAIL  UPPMSI F T ,MA X ) 

CALL  DPPHStF  P.MAX) 

CALL  MUYSRC (Ftl,FIP,CnSP,SINP) 

CALL  DPPHSt r IT ,MAX ) 

C A 1 L npPHS  ( F TP  ,M  AX  ) 

Pi  FITS  I t S ,N  I = WF  AL  (IT) 

PLiiTSI  (ft,N)  = R(  AL  (IP) 

PIOTST ( I ,N )-RF  AL ( F TT ) 

PLOT  ST ( 6 ,N ) = RF AL ( F IP ) 

If  ( THF  I A .r.T  .on.  ) WRITF  (iS,Pfl)  THFTA  ,FTT  ,ETP 
FORMAT  ( F 0.  ? , I IOC,  ?(  F8  .?  ,F  7. 1 ) ) 


CALL 

CALL 

CALL 

call 

CAIL 


rr  f . I I 
6,80) 


IF  (THFTA. I F. 00.  ) jVHIIF 
t,F  TP 

format  ( F O. ? ,8(F 8. ?,F7. 1) ) 

COOT INUF 

iMXXil  hoof  ♦ I NTHF  T A-1  ) »01HF  T A 

FAIL  PL  1 ( K I (MAX , FOR ,07, PHI ♦ I 80 ./PI  ,FRFO/IO( 

A ) 

font  iNu) 

FONT  I NUF 
GO  IF)  S 
STOP 

I on 

SliPROUIIOt  PITCX  (TMAX,FOD,niA,PHl,EREO,NTH( 
01  ME  MSI  ON  0AIA(S),JTITLF(?0),PLF)TST(P,)80) 
DATA  IFIHST/l/,NCrM/7/ 

OAlA  XAXIS/0./,  YAXIS/6./,  VL/-77./,  VT, 
C W ; 1 . ♦ HI 
VS--7.*  HI 
1 X : 10  . /0I 
Rf An(s,  ) I 
forma  I (MX 
IF  l.O  CASF 
IF (NTYPF 


■I 


S Rf  AO  ( S , 1 I NT  YPF  , OAT  A 
I F ORMA I ( MX , I 7,SF 1 0.0  I 
RF  UiRN 

.fO.  ?)  R11URN 
If  ( I F I R s T . ( 0 . F I (-0  If)  ? 

first  PlOTTlriF,  CAtl  INITIALIZF 

F A(  I PI  (IT  I 7 . , 1 .,  -I  ) 

I F I M I YPt  .N(  . I ) GO  TO  3 


OS  u 

4040 

O'lf 

40  so 

O'U 

401^-0 

ns  If 

40  ^0 

nsu 

4 0 PO 

hmj 

40QO 

osu 

4 100 

nsu 

4 no 

f,SU 

41  ?o 

f:sii 

4 1 10 

osti 

4 1 40 

nsu 

41  so 

ns  u 

4 1<S0 

nsu 

4170 

nsu 

4 1 BO 

nsu 

4 1 90 

nsu 

4?00 

nsu 

4?  10 

nsu 

4??0 

nsu 

4?10 

nsu 

4?40 

nsu 

42S0 

nsu 

4?O0 

nsu 

4?70 

f«Su 

4280 

USlf 

4?90 

nsu 

4100 

nsii 

4110 

nsu 

4120 

nsu 

4110 

csu 

4340 

nsu 

43S0 

nsu 

4100 

nsu 

4170 

nsu 

4180 

nsu 

4190 

csu 

4400 

nsu 

44  10 

nsu 

4420 

nsu 

44  30 

nsu 

4440 

,EFBT ,FEBP  ,ET 

rFP  *FT  TOSU 

44S0 

osu 

44f  0 

nsu 

4470 

nsu 

44B0 

nsLt 

44  90 

NTHFTA,PLOTST 

,nTHFTAnSU 

4SOO 

nsu 

4S  10 

fJSU 

4S20 

nsu 

4S10 

osu 

4S40 

nsu 

4SS0 

nsu 

4SO0 

PLOTS! ,DT ) 

nsu 

4S70 

nsu 

4 s H 0 

nst» 

4SO0 

, PC/0./,HT/. 

lA/  nsu 

4600 

n<u 

4610 

nsu 

46a’’0 

nsu 

4<  30 

nsu 

4^  40 

nsu 

46S0 

r>su 

4660 

nsu 

4670 

nsu 

*tOPO 

f^SU 

4/90 

nsu 

4 ?00 

nsu 

4710 
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'4 


c A«is  CASE  nsu 

IF  tIFTRST  .fO.  0.1  CALL  PLOT  I XA  » I S ♦ , 0 . . - I OSU 

If  ( OAIA(l)  ,NI . 0.1  XAX!S=nATA( I ) OSU 

IF  ( OATA(?)  .NC.  0.)  VAX  IS  = OATA( /’)  OSII 

IF  I OATA(l)  .NF.  0.)  VL^OAIAni  fISIJ 

IF  ( OAIA(A)  .NI.  0.)  VT=;|)AIA(A)  OS  U 

IF  ( oatais)  .NI . 0.)  pc^i.  nsii 

IF  ( oAiAisi  .Ft'.  0.)  PC=-o  nsu 

(,  PC  IS  RFAL  CtniF  VAHIABLF  I.  IS  INPIVIOIIAL  PHI  LABELS  OS U 

OXiTMAX/XAX I S nsu 

nv-ivT-vi  i/vAxis  nsu 

i(  pv-o  nsu 

CALL  Axisto.  ,0., ’tufta  idfc.)  • ,-i  1 ,xaxi s,o.  .0. ,nx , i . I nsu 

CALL  AXIS(0.,0.,  •F-FIFLU  ( (IB  I * , I ?,  yaxi  s,«jo  . ,vl  ,oy,  I . I nsu 

C DRAW  f/l),  niA  AND  fPFO  LABELS  OSU 

XCf'L  = XAXI  S-F  . QSU 

ycol=yaxi s-.?s  nsu 

XC  =XCUL  ♦ S . ‘LW  nsu 

call  sympul  I xc  ,yc('l  ,mt  , • F/n=  • ,0.  ,a  I nsu 

call  NLiMBiKixc*  6.  »c  w , YcnL  . MT,  Enn.o. , 3 I nsu 

YCnL=YC.nL-vs  nsu 

CALL  SYMBC’L  I xc  .YCOL.ur  , *01  A OF  REFL  ( I N I = • , 0 . , 1 7 ) nsu 

CALL  MUMBF P ( XC ♦ 1 H .»C W, YCOL  ,HT ,ni A , n . , ? I nsu 

YC  (iL  = YCnL-VS  nsu 

call  SYMBOL ( XC ,YCnc ,HF ,• FREQ  ( GHZ ) = • , 0 . , 1 1 I nsu 

call  numpfr ( xc « I ?. »cw , Ycnc ,hi , frfo ,0. , ? i nsu 

YCtiL  = YCnL-vs  nsu 

IF  ( PC  . F 0. 1 . ) GO  TO  A nsu 

(.  PR  AW  COMMON  PHI  label  nsu 

CALL  SYMBOL  I XC  ,YCriL  ,HT  ,•  PHI  CUT  U)  FG  ( = ’ . 0 • t I A I OS  U 

CALL  nlimpfr  ( xc  4 1 s. »cw ,YcnL .HT , PHI , n.  ,-n  osu 

YCOL=YCOL-VS  nsu 

GO  TO  A nsu 

T IE  1 ntypf  .Nf  . A)  c-n  in  o osu 

c PLOT  title  nsu 

NC=PATA( I ) nsu 

RE  API  5 ,7  ) I 1 1 1 1 F nsu 

7 format (70 AA)  nsu 

XT r ( X AX  I S-NC XCW)  /? . nsu 

CALL  SYMBOL  ( X I .7'.,H1  , I I I TLF  ( 1 I ,0.  ,NC  I OSU 

on  TO  A nsu 

6 IE  (NTYPF  .NF.  ?l  GO  TO  S OSU 

C PLOI  CURVE  nsu 

IF  ( OATA(I)  .NF  . O.I  NCLlMrOATAin  OSU 

ICRV=ICRV-*1  nsu 

NPrNTHETA  OSU 

IF  ( NCO‘*.LT.  S .AND.  TMAX  .GT  . 90.)  NP  = 90./0T*I  ns  11 

01)X=XAXIS/(N1MFT*-1  I nsu 

I C ^ ’ nsu 

no  fl  l-i.NP  nsu 

X- 1 i-i )Aonx  nsu 

Y-  (PLOT  SI  (NCOM  , I I- VL  I /PY  nsu 

IF  ( Y .LT.  0.1  Y-0.  nsu 

IF  ( Y .GT.  YAXIS)  Y=YAXIS  OSU 

CALL  PLOT  (X  ,Y,IC  I nsu 

IF  ( ( I/I  X )♦  IX  .F C . I I CALL  symbol ( X, Y, .07, ICRV.O. .-I  I PSU 

p ic^?  nsu 

C LAPEL  CURVE  nSU 

call  SYMPni.  ( xc -1  . , YcnL4Hi /7  . , .07  ,icrv,o.  ,-i  > nsu 

(ALL  PLOI  ( XC- .7,  YCOL  *HI/?.  . zi  OSU 

call  SYMPPL ( XC-. Z.YCni *HT/7. , .07 ,icrv.o. ,-I I nsu 

CAIL  L ABI  L ( XC  . YC  PL  ,HT  ,NCOM  , 'PR  INC  STRUT  SCAT  •,  OSU 

1 'X-POL  MIMII  SCAT  •.•PRINC  Effp  gCAT  ','X-PLU  FFEP  SCAT  OSU 

7 *pr;nc  fifl  * FFin*,  •x-pol  rfel  * effp*  , 'priaiC  total  fiflo*,  nsu 

A ’X-POL  T(>TA1  flfLP'l  P'u 

IF  ( PL  .FO.  0.1  GO  TO  9 nsu 

C LABEL  INOIVinUAt.  PHI  CUT  OSU 

fail  SYMPni  ( XC  4 1 p.  •(  w . yf  ni  ,ht  . • phi  I • ,0. , a I nsu 


A 7 70 
47X0 
4 7 AO 
4 7 so 
A 7 AO 
A770 
A7B0 
4700 

4 BOO 
4B  10 
4B70 
AP30 
ABAO 
ABSO 
AB60 
A«70 
ABBO 
AP90 
AOOO 
A910 
A970 
A930 
A9A0 
A9S0 
4960 
49  70 
A 9 BO 
4990 

sooo 

SOlO 
SC)  70 
SO’.O 
S040 
SOSO 
S0<  0 
so  70 
SOPO 
S090 
SI  00 
SI  10 
SI  70 
SI  ?0 

5 1 40 

SI  so 

SI  AO 
SI  70 
SI  BO 
5190 
S700 
S710 
S770 
S770 
S740 
A?SO 
S7A0 
S7  70 
S7B0 
5 790 
S300 
S3  10 
Sx^O 
S3?0 
S340 

S3  so 

'3  AO 
S37P 
S3P0 
S 300 
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t At  L V t XC  *7  *CW,  YU)L  . H T , PH  1 , U . , - 1 ) 

nsu 

S4.0(t 

YCOl  sYCm  -VS 

nsf/ 

s*;  fo 

( A U P L f If  1 .> 

('Sti 

sA?n 

H 

IF  1 R S T = 0 

nsu 

* AlO 

f.h  Tn  s 

fiStJ 

Si.  An 

f NO 

nsu 

SA  sr- 

SUPRIUJI  !Nt  LAPK  (X,Y,H,N,Cf  ,C2tC!t,r^,CS,C^ 

,f  7,f  R 1 

nsu 

sA^.n 

Di Mt  N‘  ION  c 1 (1 ) ,c ? 1 1 ) ,c^  ( n ,c^  ( n ,c^  I n ,C6(  n ,C7  n ) ,r.fl ( 1 1 

flSU 

SA  70 

f.n  in  ( 1 , »,H)  .N 

nsu 

SAflO 

cut  svMfioi  ( X ,y, 1 1 1 ) ,0. . 1 7) 

nst> 

*,A<iO 

&(1  TO  9 

rtS(/ 

ssno 

7 

r At  1 <;yMpaL  ( < , Y,n,c?(  1 ) ,0. , 1 7) 

nsiJ 

s*  10 

f.n  in  >7 

nsw 

C M 1 SYMHPl  ( X . y,  M,C  1 ( II  ,0.  , 1 7 1 

ns  u 

SS'iO 

cn  TO  *' 

nsu 

SSAO 

<• 

Cau  syMpnL(x,y,H,rA(i(,o.,4  7) 

ns  11 

ssso 

&fi  in  9 

nsu 

SS7.0 

Cau  SYMBoi  ( X,  v,H,r 1 1 ,0 . , rn 

nsu 

SS70 

GO  in  V 

nsu 

ssno 

C A Li  '.YMiuii  ( X , y , H,r. M 1 ) .0. , 17 ) 

fistf 

SS  VO 

ffi  rn  9 

nsu 

• ^.on 

CAU  SYMPflL  1 X ,y,H.C7(  n ,n.  , 17) 

nsu 

S<S  10 

fn  in  4 

nsu 

so?o 

CAU  ^YMpnu  X , y ,H  .c  R n 1 ,0. , 1 7 1 

nsu 

‘A30 

<1 

Kt  Turn 

nsu 

S7.A0 

1 Nl) 

osu 

ssso 

MlflRUllTINF  Ml  lYSBC  11  1Ht  1A,FPHI,CnSP,SINP  I 

nsu 

ns(i 

SF.  70 

PnuiINr  CONVIRT  ITHIIA.CPHI  to  f huygfn 

PRINCIPAL  AND  CROSS 

nsu 

s/,ho 

P(UAR17Arn  COMI’ODFNTS 

nsu 

Sf.VO 

ns(/ 

*>700 

riiMi'i  F X f ThF  I a ,F  Fmi  , ‘KIP 

r^st) 

S710 

STnRiCflSP»FTHf  lA-SIMP'fPHi 

nsu 

S 7?0 

F 1 HI  1 A - S 1 NP  •!  I HF  TA  «f  OSPAFPMI 

nsu 

s ryo 

I PHT-bKlP 

n<:u 

S7A0 

Rf  TURN 

nsu 

S7‘  0 

FIJI) 

nsu 

S 7^.0 

SlIURinn  INF  ‘1RUT(PnI,IHFlA,PHlS,7?,R5i’,RLn 

,rsT,FSPi 

f>SU 

*>7  70 

nsu 

STi^O 

ns(/ 

S700 

DiBiFN' ION  Rf,y(ino),ASV(iooi.R;yMoo),A7yiino),iFMOo),nFHoo) 

nsu 

SAOO 

ClIRPLFX  J,C.0,C()1  .1  0?,rPR.ns,l)H,T|MP,1f  RMS, 

TFRf-tH,  AMf  ?S  ) ,0Y, 

nsu 

*yf\\0 

FSllMY,  t AT  ,FAP  ,FOT  ,FGP,EALP,FGAM,CP(),FST  ,FSP 

,OYn 

nsu 

Sfl?0 

COMMON  /SIR/J.GPn.AM.nM.PI.RK.NX.nFL.A.TF, 

OF  ,BFT 

nsu 

SR  30 

IMl  =7HF 1 A*P I /I RO  . 

nsu 

ShAO 

cii.'.r^cnsi  iHi  ) 

nsu 

SftSO 

:■  INT  = SIN  1 Tt'F  I 

nsi; 

iKXPS-SKKPNJS) 

n<  ij 

^ R70 

CFi5p‘;  = cns  ( PHI  s 1 

nsu 

'PRO 

S FNR=  INI  111  It 

nsu 

*<R^0 

ct)'.R‘rin,(M  1 1 

nsu 

*^000 

f niRrCDSH/SK.'P 

nsu 

SVIO 

?nr7?«RS?*rmh 

r«su 

RSI-RS?-PUl»GINR 

nsu 

SV30 

IF  I«S?.(.T.At  PC7=A 

rsii 

S<^AO 

PU'IAI/SINR 

nsu 

* ‘iSO 

RP^RSP/SINf) 

nsu 

L 1 = RS  I /nu  ‘ ? 

US(i 

70 

ty=P‘?/OI L‘ 1 

nsu 

*=  V RO 

N=l  ?-l  1 

ns(j 

NI  -N‘  I 

nsu 

7.noo 

ni  1 S - 1)1  1 / SI  NR 

nsu 

7.010 

RL  1 - 111  L S Ml  1 - 1 I 

n«-  u 

7.0  ?0 

Ri  ?^nF  1.  s*  1 1 ?-n 

fiVU 

7-030 

CP'  iC  F xp  1 J*  1 p I /R  .‘RF  ♦ ?n»  1 cnsi-1 . 1 1 1 

rSii 

7 fiAO 

71  ( XP  1 .1*11  / A . I/I  PI  *SINrt) 

USu 

7 C"0 

PH  1 ii.  PHI  Phis 

f*r.  u 

7.0*  0 

cnsfiK  ns  1 PI  1 1'l 

nsu 

on  rn 
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10 


1? 


in 


10 


?o 


S IN'I^S  I'J  ( PMI  11) 

c ns* r SINT*?  iNf*«r. nsn-cnsT*cnsR 

*l  P*ARC(t*>  IC(I*’»  1 
S I NA  = S IN  ( AL  I’  1 
>;sf.=  s I NT  * SI  Nh 

sr  r,T-MN  1 ♦(.(>?  p *c  nsn-cnsT  *s  INB 

(.AM=  *7  AN?  ( SSl,  , sc  (■  I 
S I NG  - S IN  I I.AN  I 
CrSf.=  C()S  I GAN  ) 

ALPHA  - AL  P*  I PO.  /p  1 
r.AMNA^LAN*!  nn./p  i 
0S  = (O. ,0.  I 
DM- < 0 . ,0  . I 
n:i  10  Mil  ,Mx 

CnSN-Cfl'  I M - I)  ♦ ( PI  -GAM  I 1 
7 F SNS=  2.  • AN  ( M I *C  tISM 
7 F pNMi  2 . ‘PN  ( N)  *C.  ns“ 

ns  = ns«  7F  F!MS 

nWrOH*  7FPNM 
CONI INUF 
p>.iUNP»ns 
0H=  7 f NP*  OH 
( Ap-cnsA •COSH 
RL=«L I 

no  1?  1=1 ,Ni 

CO=CfXPI.IN«K*Rl*CAHl 
I*L 1-1 

F I =-(  1 F ( K i-or  ( K I •c.nSPSi 
siiNY=F  i*r,t' 

RSY ( 1 1 -R7  AL ( SUMV  1 

A S y ( I I r A I NA(,  ( SONY  » 

RL=RL»nF  LS 
CONI INOF 

CAIL  OSF (OF LS,RSy,R?Y,NI I 
call  OSF (pr LS,ASY,A/Y,N1 I 
CYn=-R;Y(Nl  )•  J*A/Y(N1  1 
coi=cexp(j*PK*(pi«RLii*CAP/?.i 
C0?=CFXP(J*n‘*(R?«RI?l*CAP?7.| 

FNl=-IPLl-''ll*MFILl)«7Mll-ll-inF(Lll*0F(LI-I)ll/r, 
fN?  = -(R?-RL?)*(TF(l?l-F)F(L7ll/7. 

uY  = 0Y()4r  N 1 »coi  «F  N7 »r.o? 
F*LP=-CP‘^*nS*SINA*OY*SINPS*CPO 
FGA>A=CPS*OH*SI  NA  *ttY*  ( -COSOS  1 »CPO 


nsu 
ns  I 
ns  I 
nsi 

nsi 

ns,  I 

nsi 

nsi 

PS  u 
F’M 
ns  I 
nsi 
ns  I 
nsi 
psi 
ns  I 
nsi 
osi 
ns  7 
psi 

ns  I 

PSI 

nsi 

nSL 

PSI 

nsF 

PSI 

nst 

OSI 

nsi 

nsi 

osu 

nsi 

nsi 

nsi 

OSI 

nsi 

PSI 

nsi 

OSI 

nsi 

nsi 

nsi 


F A I = -( roS**CnsG*  CnSB«SINA*S INB 1 *CnS 7»Cnsn* COSA * S ING»CnST »S I NO* S I N70S L 
7*(CnSA*CnSf.*SINP-SINA*C0SP  1 nSL 

F AP=( siNA*s iNp.cnsA»cnsG*cnsp I *s iNn»coSA*siNG*coso  nsi 

FGI  = SING*  If  nSP*COS7*CnSO-SINP*SIN7  1 <Cn SG*C0S7 * S I NO  OSI 


FGP=cnsG*cnsn-si ng*c  nsn*siNO 

F A I = r ALP*  F A 7 
FAP  = (ALP*FAP 
FGP-  FGAN4FC.P 
FGT=FGAN*fG7 
F S7  = F A I . F GI 
F SPrF  AF'4F  GP 
RF  TURN 
FNP 

Sl)P«nU7  IMF  LASOIPSI.PX.FF.M 
ni  «F  NS  ION  PS  I ( so  1 ,PX  ( so  1 , F F ( so  I 

WR  II  r (6,7) 

FOPNAT  (//MG, '0  IFCFWISF  LIN'AR  SFIP  I MPI|7  • , //I  1 P , ’PSI  • , 7 SI 
on  70  I - 1 ,N 
RF  A0(  S , I P 1 P>  ( I 1 ,F  F ( I I 
FOPNA  7 ( I I 0. 7 ,r  I O.  s ) 

FFIRN*  I ( 7 U',F  10.  ,',f  IS.SI 


WR  I 7f  ( / , 
IF  ( I .(  C. 
P S I ( I - 1 I ■ 
CON7 INOF 


!•  I PX(  I I ,FF  ( I I 

1 1 r.n  III  70 

( I F I I I- r F ( l-i 1 I /( PX( I ) -py ( I-I n 


nsL 
ns  I 
nsi 
nsi 
ns  I 
nsi 
nsi 
nsi 
nsi 
nsi 

nsi 

nsi 
' ,/insi 
nsi 

OSI 

PSI 

rsi 

P'l 

PSI 

nsi 

nsi 


I 

i 

f 


<K0«0 
^>090 
MOO 
f \ 10 
ftl?0 
M 

^ 1 ^.O 
M ^>0 
6 W ^ 
( \ 70 
M PO 
61  90 
6?00 

6?20 
6?30 
h?i*0 
6?S0 
6260 
6270 
6?rn 
6290 
6300 
6310 
6320 
6330 
6340 
6350 
6360 
6370 
^3  >10 
6390 
^400 

64  10 

6420 
64  30 
6440 

64 

6460 

6470 

6400 

6490 

6600 

65  10 
66;>0 
66?0 
6540 
66S0 
6»‘60 

65  70 
6500 
( *'90 
6600 

66  10 
6620 
66-»0 
6^  40 
< ( •'0 
6660 

66  70 
66^0 
6 690 

6 7 00 

67  10 

7 7?0 
67  >0 
6 740 
7.  7'  0 
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oo»'^<^-‘^or'or^ooo  00*^0  f^.'^r'ooooooor'ooop.r'oof^or'r'oo 


\0 


10 

?c 


Rf lURN 
tNO 

Sim^nUT  INf  L SUM  I PSl,PX,rF.r>S!tN3,FF) 

01  Mf  N$  ION  OS  n ‘>0  ) ,PX  I *>0  1 ,F  F r>0  1 
0(j  10  1 = 1 |N3 

IF  IPSr.GT.PXt  IMM  GO  TO  10 

ff^vfii)4[)Smi  »*<pn-pxiin 

«F  TURN 
CONT 1 NUF 
f f =n. 

RF  TURN 
FNO 

SUPROUTINF  ORPMS (r ,MAX 1 
COMPlf.  X i 
RF AL  OP  ,F ASF  »MAX 
P|=3.1AI59?6*»3 
AF=CABS(F  ) 

IF  { AF  .G1  .0.1  GO  TO  10 
0P=-*‘00. 

F A^f  =0. 

GO  TO  ?0 

0P  = 20.AALf)G10(  AF  1“MaX 

FASFnAlAN?(AlMAGlFl»RFAl.(F))AlflO./PT 
F=CMPLX(OP,f ASF ) 

R*"  UIRN 
FNO 

SUnPOUTlNF  GTOIPHI  ,PHIP,PFTA0,S  , S P , WA  , S L * A S , I LOP  * OS  * OH | 

PIJRPOSF  - TO  COMPUTF  OlFFPACTION  COFFFlCIFMTS  OS, OH 
PAPAMETFRS 

♦ ♦A  inputs** 


PHI  S OIFFRACTFO  ANGIF 

PM  IP  ^ TNfrf»ffir  anglf  at  thf  fogf 

BFTA0=  ANGLF  PFTWEFN  THf  TNCTDCNT  RAY  THF  TANGFNT  70  THF  FOGF 
S = OlSTANCf  from  FTFLO  POINT  TP  THF  FOGF  ALONG  THF 

otffractfo  ray  path 

SP  = oistancf  from  Thf  sourcf  point  to  thf  fogf  along  thf 

JNCJOFNT  ray 

WA  S WFOGF  ANGLF 

SL  = oistancf  PAWAMITFR  WHFN  FXTFRNALIY  PROVlOFO  (WHFN 
ILPP  = CflNTROL  fniMfH  p 


I LflPr  1 
UOP  = ? 
UOP=  3 

ILPP-A 

• ♦ AOHT  PUT ♦ 


PLANF-WAVF  INCIOFNCF 
: CYLINPRICAL  WAVF  INCIOFNCF 
: SPHIRICAL  0»  CONlCAL-WAvr  TNCIOFFiCF 

thf  oistancf  paRAMFTfR  L-SL  TS  pROVTOFO 


OS  I 
OS  II 
ns 
ns 
ns 
nsu 
ns 
os 

PS 

ns 
ns 
os 
us 
ns 
ns 
ns 
ns 
ns 
ns 
ns 
ns 
ns 
ns 
ns 
ns 
ns 
os 
ns 
ns 
ns 
nsv 
os 

ns 
ns 
ns 
ns 
os 
os 

r‘SL 

os 
ns 
os 

PS 

TLrP=AOS 

os 


AS 


OS 

OH 


s AMPlTTlIOF  OF  THF  FIFLO  AIONG  OirrPACTFn  RAY  (NOT 

WHEN  UPPs^l 
COFFF  JC  IFNT 
COFFF IC  n NT 


= SOT  T 
r HAR[) 


DIFFRACTION 
OIFf  P act  K'N 


ns 

ns 

os 

FXTFRNALLY  f'S 
PS 

ns 

ns 

APPL ICAOS 

ns 


♦••UNITS*AA 


ANGLF  t nroMif 

LINGTh  : wAVILINGTh 

RfOHlPFO  niN( T IONS  t TRANSITION 


FliNfTInN  F(X1  ANO  SIGN  FUNCTION 


f> 

r*S 

ns 

ns 

os 

ns 

ns 

SGOS 

r^S 


RMFRCrjCI:  A IlNirOPM  CfOMfTPICAL  THFOPy  nr  nirfRACMON  FOR  AN  rncrosi 
IN  A PfRMCTLY  CriNOUCTING  SURFACF  , PY  Q .G  . OHYOl  m J I AN  ANO  P.H.  OM 

paimax  , ppnr.  of  iiff  , vni  nov.  irta  r»si 

ns  I 

ns* 


A760 
6770 
f 7P0 
f 7<^0 
6M00 
6R10 
6R?0 
6R30 
6RA0 
6PS0 
6n60 
6fl  TO 
6R80 
6R90 
f QOO 
6P10 
6S»?0 
6930 
/9A0 
6OS0 
6960 
6970 
69P0 
6990 
7000 
7010 
70P0 
7030 
70A0 
7OS0 
7060 
7070 

70PO 
7QO0 
7100 
71  10 
7120 
7no 
71  AO 
7 I SO 
7160 
7170 

71  PO 
n 90 

72  00 
T?10 
7720 
7230 
7?A0 
7?S0 
7260 
72  TO 

72  «0 
7290 
7300 
7310 
7S?0 
7330 
73A0 
7360 
7360 

73  70 
7*^80 
7 390 

taco 

7A  10 
7A  20 
fi,  AO 
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r oMPi  f X t , A .1  ,co,  1 1 ,c ? ,c  t ,r  A ,f)S , DH,f  4S^  ,0 

I,-*  h 

Rt  AL  L ,K  1 ,K2 

r/  u 

7A*.r 

A(l,WN,N)=/.0»fl)5.(».lAl‘i<J?<,S»WN»N-'T/?.(1|**? 

n^u 

7a  ( 0 

( 

osv) 

fu  70 

c 

* « • 

CnwPUIF  rjPrpSSARY  CIINSfANTS 

nsu 

74  MO 

r. 

OSil 

74  VO 

PI  = 3 . 1 A 1 

nsir 

7S00 

P M=^.0»PI 

TS\0 

PlArPI/A .0 

osu 

7sro 

PI  1RO  = PI/Ift0. 

r»su 

7‘‘>30 

A J=CMPl  X ( 0.  , I . ) 

nsu 

7S40 

t PSN  .ooni 

f'SU 

7SS0 

C'SU 

?snO 

GA=1  . 

rsu 

7S70 

GH=  1 . 

nsu 

7880 

WG  = 3pO,~WA 

nsu 

7SV0 

c 

nsu 

77>00 

IF  ( PM  IP  . t 0 .0.  .flH  . PH  IP  .f  0 . WG)  GS  = 0. 

nsu 

7010 

IF  ( PHI P. LU.O. .OF. PHIP.FO.WG ) GH=.5 

nsii 

70?0 

c 

OSH 

77a  ■^0 

c 

* • * 

COMPOTF  ANGIFS  RlTA-fFRn.PFTA-PO^ITIVE.BFTA-NFGATIVF 

nsu 

7f>40 

c 

nsu 

77>SO 

WN-W('./  HO. 

nsu 

7600 

hf  TApsiPHl4pHlP|*PllH0 

osu 

76  70 

8f  T AM= ( PM I-PHI P) ♦PI  1 80 

nsu 

7/.H0 

c 

nsu 

76*^0 

c 

• • « 

CHOOGF  DISIAfICf.  parameter  FOR  TYPE  OF  EGOE  ILLUMINATION 

osu 

7700 

c 

nsu 

7710 

IF ( ILOP.TO. 1 1 L=S*SIN(BETAO»PI 

nsu 

7 7?0 

IF  ( ILOP.FO.;>)  1 = 5,»SP/ ( T,*SP  ) 

nsu 

7730 

IF(IL0P.fQ.3l  L=^S*SP/(S«SP|*SINIRFTAO*PII80)»*? 

osu 

7740 

IF  ( IlOP.f  O.A  I 1 = SL 

nsu 

7750 

c 

osu 

7760 

r 

* • « 

COMPIIIF  A(r.,LPI  IMF  AMPLIUIOE  OF  THE  FIELD  ALONG 

THr  OlFFRACTEt) 

RAOSU 

7770 

c 

nsu 

7700 

IMIlOP.FO.il  ASrl.O/SORTIS) 

nsu 

77«>0 

I F I ILOP.FO.?  > A5-I -O/SORT ( S*SIN (BFTaOaPI IRON 

nsu 

7fl00 

IF(IIOP.EO.T)  as=S0RT(SP/IS»ISaSP1 1) 

osu 

7M  10 

IF  1 ILOP.FQ.A  ) AS  = 1 .0 

nsu 

7H?0 

c 

nsu 

7 130 

c 

« • • 

COMPUIl  ANGUl AR  arguments  FOR  COTANGENT  FUNCTION 

nsu 

7P40 

r 

nsu 

7fl60 

ANG1= (PI  ♦flFTANl/12.0*WN) 

fisi; 

7060 

ANG?=(PI-R( TAN1/(?.0»WNI 

nsu 

7H70 

ANG3= ( P I « BF  I AP 1/ (? .0*WN 1 

osu 

7080 

ANGA=(PI-HF  lAP  1/  (2.0»WN) 

os  u 

7PV0 

c 

ns  VI 

7V00 

r 

• # « 

calculate  N«  ANO  N- 

nsu 

7V  10 

r. 

osu 

7020 

NPN=  IF  IX  ( (PI  -.BF  T AN  1/  ( P I I *WN  1 *0.  A 1 

nsu 

7V30 

NNN-I F IX(  (-PI«Rf  TANI/(PII*WN(«O.A) 

nsu 

7V40 

NPP=IFIX((PI*RFTAP)/(PII»WN)«0.A) 

nsu 

7VS0 

c 

nsu 

7w^.O 

c 

« • • 

COMPUIl  API.UMENTS  OF  TWANSITION  FUNCTION  FIX  I Xs 

Kl  A 

nsu 

Vt  70 

ARGM  P I I AL*  A ( BF  T AN  ,WN  ,NPN  1 

usu 

yohO 

AKG?  = P11»1*A|RIT  AN  ,MN  ,NNN  1 

0'  u 

7V40 

A0G3-PII*L»A(RFIAP,WN,NPP) 

nsu 

Mono 

ARGA=P I I AL ♦ A ( BF I AP , WN.NNP 1 

nsu 

MO  1 0 

r 

nsu 

MO.'O 

c 

* ♦ • 

(iS  ANO  OH  CAiCULATin  BY  XFUFR'S  FORM  WHEN  FIFLO 

pnun  IS  NOT  IN 

ns\» 

8030 

NNP  = I F I X I ( - P I . F(F  1 AP  1 / 1 P I 1 »WN)»  0.  S I 

usu 

8040 

r. 

transition  RFGIONS 

nsu 

80  SO 

c 

(>SU 

HO^  0 

IF  { AHGl  .Lf  .W)  C,r\  If)  *)0 

nsu 

Mr  ? - 

IF  ( Awr,?  .1  f .W)  GfJ  TP  ^0 

u 

8r^r 

n ( APr.^.  1 1 . w ) Gf‘  in  so 

•'*  II 

M0'*«- 

1 f ( AP(,S  . t f .W  1 f.O  Tf»  so 

nsu 

M 1 ' 

rOiSIMP|/wNI»(fXPI-A.)«PIAI/(WFJ*PII»SIN(R(TA0*PIlH0ll 

T 1 t 

^ M 

XO-A043  400 

UNCLASSIFIED 

20=2 

■•••oo 


OHIO  STATE  UNIV  COLUMBUS  DEPT  OF  ELECTFICAL  ENOINEERINS  F/0  0/1 
*iif^2^*^***^***  CIRCULAR  reflector  ANTENNA  INCLUDINO  FEED— ETC (U) 

AU«  77  s M LEE#  « C FUOOUCK.  C A KLEIN  F30602-76-C-02a4 


4SB1-1 


FAOC-TR-77-tSO 


9 -77 


L 


r*  o r»  o o 


c 

c 

c 

c 


c 


c 


c 


c 


c 

c 

c 


c 


K 1 = 1 .(>/( rns t PI /wN(-cos( AN/wNi  i 
K?  = i .o/(c  n«,  ( p|/wNi-cns(  Rf  TAP/wN)  i 
i)S=f.5.*tn»  IK  i-K?) 

DH=GH*rri*  (K 1 <K  ;>) 


RF  TURN 


CAICUIAII  THf  FOUR  TfHMS  OF  THF  OTFFRACTTON  CDFFFlCin 


SO  Cn=-CFXPI-AJ»PIA)/(?.o»WN*Pn*SIN(RfTAO*PIIflO)l 
IF  IABS(ANf.l-Pl*NPM).GT.FPSI  ) GO  TO  100 
roS*l  ANi;i-PI»NPN  )»2.»WN 
C1=0IFPS,WN,L.P1  ) 

CO  TO  ISO 

100  r,  1 if riTANi ANCi  1 »F  (AR&n 

ISO  IF  (AHS|ANf.2*PI*NNN).GT.FPSn  GO  TO  200 
FPS=(ANG2«PI AUNN )*?.»WN 
C7=0lf PS,«N,L,PI  I 
GO  TO  2S0 

20F  C2=CnTAN(AMG2l*F (APG?) 

2S0  IF  I AP  S I ANC.  T-PI*  NPPI  .GT  .FPS  1 1 GO  TO  300 
L PS=I ANG3-PI »NPP  1»2 .AWN 
C3iOltPS,WN,l ,PI  I 
GO  TO  3S0 


300  C3=C0TANI ANG3IAF (ARG31 

3S0  IF  I ABS I ANGS aPIANNP I ,G I .F PSI I GO  TO  AOO 
F PS=I ANGSaPI anNP  1 a 2. awn 
CS=OI f PS.WN.l ,PI  I 
GO  TO  ASO 

AOC  CAifOTAN I ANGA I *F  lARGA  I 

ASO  continuf 

aaa  computi  os  ano  dm 


1 

2 
3 


nSif.SACOA  If  lAf  2-C3-CA1 
OHiGHACOAlC I AC2AC3AfA I 
RF  TURN 
FNO 

function  SGNIXl 

PURPObf  : TO  GIVF  MINUS  SIGN  IF  X IS  NfGATlVF  OR  VICF  VERSA 
IF  IXI  1,2,3 
SGN=-1  . 

RF  TURN 
SGN=0. 

RFTURn 
SGN= 1 . 

RFTURN 
I NO 

COMPLEX  I UNCTION  F IX  I 

transition  function  FIXl  FOR  DIFFRACTION  COFFFICTFNT 


FIX)  IS  COMPUTFD  BY  FRFSNFL  INTFGRAL 

FOR  X>h.,OR  XI.7,  FIX)  IS  FVALUATFO  BY  APPROXIMATION 
FRROR  BY  APPROXIMATION  IS  WllHIN  II 
SUBROUTINf  REOUIRCO  CSXIC,S,XI 


COMPLFX  IX, t 

RFAL  PI/3.IAI5‘»2ftS3S/,PIA/.7BS39Rlt/,  PQ  / 2.S0A>S?B3/ 
COMPlFX  CPI  A/ I . 7 07 1067 R, .707 1067 fl|/ 

IF  IX-P.I  20,20,10 

in  x?-x*x 

X 3 = X?AX 
Xa=X2»X2 


nsu 

iisu 

os  II 


nsu 

osu 

osu 

flSU 

osu 

osu 

nsu 

osu 

osu 

nsu 

nsu 

osu 

nsu 

osu 

osu 

osu 

osu 

nsu 

nsu 

nr.u 

osu 

osu 

nsu 

nsu 

nsu 

nsu 

nsu 

nsu 

nsu 

osu 

nsu 

nsu 

nsu 

nsu 

osu 

nsu 

nsu 

osu 

nsu 

nsu 

nsu 

nsu 

nsu 

rsu 

osu 

nsu 

nsu 

nsu 

osu 

O' II 

osu 

osu 

osu 

osu 

osu 

osu 

osu 

osu 

osu 

nsu 

osu 

osu 

osu 


91 


BI20 

61  30 
B1  AO 
RISC 
RI60 
61  10 
PI  RO 
0190 
0200 
R2I0 
0220 

62  30 
P7A0 
02SO 
6260 
R270 
R2H0 
R290 
R300 
0310 
R320 
R330 
63AO 
R3S0 
R360 
6370 
0300 
R390 
6A00 
BAIO 
flA?0 
0A3O 

RaaO 
OASO 
RA60 
RA70 
RA80 
6AP0 
0600 
R6  10 
R6?o 
flS  30 
OSAO 
R6S0 
0660 
6670 
6600 
06<i0 
6600 
0610 
06  20 
66  30 
06A0 
0660 
0660 
0670 
R6II0 
66  RO 
R700 
0710 
0770 
0 7 30 
07A0 
6760 
6760 
0770 
0 7R0 
8790 


o r>  n rt  r>  oo  o 


f =t*»Pl  X 1 1 . /X,’*'.  1.  7^>/Xii,  .S/X 

flSU 

8800 

RUDRN 

ftSIJ 

88)0 

IF  (X-.?  t 40, 3C.  10 

nsu 

8870 

ZX»CMPIX«0. ,X> 

nsu 

88  30 

C»LL  rsx(c,s,x» 

nsu 

8840 

C 

nsu 

8880 

Sl=.S-S 

nsu 

8860 

F = iORT  ( X MCI  XP(  /X  MP0»CMPLXISI  ,CI  1 

nsu 

8870 

RF  TURN 

nsu 

8880 

Z-SQKT(PI*X)-7.*  X*C P 14-7. »X*X/CP 14/3. 

osu 

8890 

f =/*Cf XPICMPIX 10., Pl4« X 1 » 

nsu 

8900 

R|  TURN 

nsu 

8910 

tNO 

nsu 

8970 

SlIBROUTINF  CSXIC  ,S,XI 

nsu 

89  30 

nsu 

8940 

PlIRPOSf 

nsu 

8950 

COMHITFS  Tuf  FRFSNFL  INTFGRALS. 

nsu 

8960 

nsu 

8970 

ofscriptiun  of  PARAMFTFRS 

nsu 

8980 

c - Tut  rfsultant  output  of  CIX) 

nsu 

8990 

<;  - THE  rfsultant  CAITPUT  of  SIX) 

nsu 

9000 

X - ARCUMfNT 

osu 

9010 

osu 

90  70 

FROM  SUDROUTINF  CSIC,S,X)  OF  SSP 

nsu 

9030 

F«AHS( X I 

nsu 

9040 

IF  (Z-4.  ) 3,3,4 

nsu 

9050 

nsu 

9060 

IF  IZ.LF.4.) 

nsu 

9070 

osu 

9080 

1 C = SORT  IZ  I 

osu 

9090 

S = Z*C 

nsu 

9100 

I- 14. -Z I * I4.«Z  ) 

osu 

9110 

C-C»M  1 M IS.1007P5F-n*Z«5.?44?«»7F-'»)*Z  + 5.451  I87F- 

7)4Z 

nsu 

9170 

l♦3.Z7330«F-5l•Z♦  1.0204  I 8F -3MZ*  1 .107544E-?  MZ*  I .8409F.se- I ) 

nsu 

9130 

S = S*(  ( I 1 16.67768  I E-I0*Z«S.R83158E-«I*Z»S.0SI141E-6MZ 

osu 

9140 

1*2  .4418T6F-4  )«7*6.1213?0E-3)*Z*B  .026490F-7  ) 

nsu 

91  50 

RF  TURN 

nsu 

9160 

nsu 

9170 

IF  IZ.GT.4.1 

osu 

91  80 

nsu 

9190 

n=fOSIZI 

nsu 

9200 

f =SIN(ZI 

nsu 

9710 

7^4. /Z 

nsu 

9770 

A=.(l  M I 1 ( 8.  768  7S  8F-4*Z-4.  I69789F -3  I *Z*  7. 97094  3F -3  >♦7-6.7978  018-3  1 

nsu 

9730 

MZ-3.09S341E-4 )♦Z♦S.9721 51E-3I4Z-1 .606478F-S)4Z-7. 

493377E-2)4Z 

nsu 

9740 

7-4 .444091E-9 

osu 

9750 

8r 1 1 ( ( 1 1 -6.633976F-44Z43.40I409E-3MZ-7.?71690E-3|4Z»7.4?8746E-3 I 

nsu 

9760 

MZ-4.07714SF-4  )#Z-9.  314910E-3MZ-1 .707998E-6)AZ4l. 

99471 lE-1 

nsu 

9770 

7^ SORT  17  ) 

nsu 

9780 

C*.S*Z*(n4*«S*B) 

nsu 

9?90 

S-  .5«74  I S4A-r)4BI 

nsu 

9300 

RF  TURN 

nsu 

9310 

FNO 

nsu 

9370 

r.flMPLFX  FUNCTUN  0(FPS,WN,RU.PI  I 

fiSu 

9330 

real  L 

nsu 

9340 

1 =RL 

nsu 

9350 

0=WN*CFXP(CMPLXIO.,PI/4.))4| 7 .4P  MSORT It lASONI EPS  1 

-4.4PIPL*FPS 

nsu 

9360 

»»CI XP(CMPLX(0.,PI/4. ) ! I 

nsu 

9370 

RETURN 

nsu 

9380 

FND 

nsu 

9390 

FUNCl ION  COTANIX 1 

nsu 

9400 

S-SINtX ) 

nsu 

94  10 

IF  (S  .NF.  0.)  GO  TO  1 

nsu 

9470 

CUIAN=1.FS0 

nsu 

9430 

RF  TURN 

nsu 

9440 

cniAN-cosix I/s 

nsu 

94  50 

RF  TURN 

nsu 

94  60 

F NO 

nsu 

94  70 

92 


MlBRniHtNt  nt  (ISII 

nsii 

(I<iu 

OS  1 1 

nsi) 

piiRPnr.F  nsu 

COHPim  IMF  J BfSStL  fUMCTION  Fn«t  * f.IVFN  ARG(I«1FNT  ANO  nBOfROSlI 

nsii 

USAGf  nu) 

CALL  PI  SJ(  X.N.HJ.O.IFR  » OSH 

nsi) 

OF  SCRIPT  ION  OF  PARAMFTFRS  nSII 

X -IMF  ARf.DMFNT  OF  IMF  J BtSSEL  FUNCTION  OFSIRFO  OSO 

N -thf  orofr  of  thf  j bfssfl  function  OESTRFO  OSU 

BJ  -TUF  RFSULTANT  J BFSSFL  FUNCTION  OSU 

0 -HFOUIREn  AfCURACr  OSU 

IF R-M SUIT  ANT  FRROR  COOf  WHFRF  OSU 

IFR  = 0 NCI  FRROR  OSU 

IFR=1  N IS  NEGATIVE  OSU 

IFR-7  X IS  NEGATIVE  OR  /FRO  OSU 

IFR=3  RIOUIRFD  ACCURACY  NOT  OBTAINCO  OSU 

IFR=4  range  of  N COMPARFO  TO  X NOT  CORRECT  I SE E RFMARKSIOSU 

OSU 

REMARKS  OSU 

N MUST  BE  GPFATER  THAN  OR  EQUAL  TO  ZERO,  BUT  IT  MUST  BE  OSU 

LESS  than  OSU 

70«10«X-X*»  7/1  FOR  X LESS  THAN  OR  FOUAL  TO  IS  OSU 

«>0«X/2  FOR  X GRFATFR  THAN  15  OSU 

OSU 

SUBROUTINES  A NO  FUNCTION  SUBPROGRAMS  REQUIRED  OSU 

NONE  OSU 

OSU 

METHOD  OSU 

RECURRENCE  RELATION  TFCHNIQUF  DFSCRIBED  BY  H.  GOLOSTFIN  AND  OSU 

R.M.  thaler, 'RECURRENCE  TECHNIQUES  FOR  THE  CALCULATION  OF  OSU 

BFSSFL  FUNCTIONS', M.T. A. C. ,V. 13, PP.102-I0B  AND  l.A.  STFGUN  OSU 

AND  M.  ABRAMOHITZ, 'GFNFRATfON  OF  BFSSFL  FUNCTIONS  ON  HIGH  OSU 

SPEED  COMPUTERS', M.T.A.C.,V. 11, l'»S7,PP.?55-?S7  OSU 

OSU 


<><.  BO 
NA<)0 
FISOO 
NS  10 
VS?0 
N530 
<>‘•40 
<>550 
<»5fc0 
V570 
<>5't0 
4S<JG 
<>F.00 
■>/.  1 0 
<>670 
<>6  10 
<>640 
0650 
<>660 

96  to 
96B0 
9690 
9700 

97  10 
9770 
9730 
9740 
9750 
9760 
9770 
97B0 
9790 
9B00 
OHIO 
9B70 
9B30 
9B40 
9B50 
9B60 


SUBROUTINE  BF S Jl X ,N, B J , 0 , I FR 
BJ=.0 

IF  INI  10, 70,20 
10  IFR=l 
RF  turn 

20  IFIX»30,30,3I 

30  IER»7 
RETURN 

31  IFCX-15. >32,32,34 

37  NTFST«7C.«I0.*X-X»*  2/3 
GO  TO  36 

34  NIFST'90.4X/7. 

36  IF  IN-NTEST  >40,38  , 3fl 
3fl  IFR«4 
RETURN 
40  IFR=0 
N I 4 N ♦ 1 
PPRFV=  .0 

COMPUTE  starting  VALUE  OF  M 

IF  (X-5.  ><>0,60,60 
50  MA=X«6. 

GO  TO  70 

60  MAi  1 .4*X  ‘C  O.  /X 
70  MB=:N4|FIX(X  >/4,7 


OSU  9P70 
OSU  9BH0 
nSU  9B90 
OSU  9900 
OSU  99)0 
nSU  9970 
OSU  99  30 
OSU  9940 
OSU  9060 
OSU  9960 
nSU  9970 
OSU  9980 
OSU  9990 
OSU 10000 
nsuiooio 
nsui 0070 
OSU 100 30 
OSUl 0040 
OSUI0O50 
OSUl 0060 
nsui  0070 
OSU10080 
(ISU10090 
OSUIOIOO 
OSUl  Olio 
OSUI0I70 
OSU  101  30 

nsuioi 40 

05010150 
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H/r8ruMAxomA  ,mr  i 

0SUIDI60 

r. 

nsuioi 70 

c 

SfT  UPPfR  IIMM  OP  M 

OSUI 01 80 

c 

OSUI0I90 

MMAXrNTf  ST 

OSUI 0?00 

100 

DO  190  H-M?f SO.MMAX.l 

OSUI0210 

f 

OSUI 0??0 

r 

SI  1 F (M)  ,F(M-1  1 

OSUI 0730 

r. 

OSUl 07A0 

FH1=1 .OF-?H 

OSUI 0750 

FM  = .0 

0SU107E.0 

ALPHA* .0 

OSUI 0770 

I K M-  ( M/  ? 1 * ? n ?0 , 1 1 0 . 1 ?0 

OSUI  07H0 

no 

.n=-i 

nsuio?<#o 

CO  TO  no 

OSUI0300 

1 ?o 

jm 

OSU103IO 

1 10 

OSU1O370 

DO  160  K=I,M? 

OSUI  03  30 

MK*M-K 

OSUI  o:»AO 

HMKs?.»fLnAT(MKl»FMl/X-FM 

nsuio'^so 

FMiFMl 

OSUI0360 

FM|=RHK 

OSUI 0370 

IF (MK-N-l  1150, 16  0, 150 

OSU103B0 

140 

8J*BHK 

OSUI  03  90 

150 

JT=-JT 

nsuioAoo 

S=1«JT 

OSU10610 

160 

AL  PHA*AL  PHA4BMK#  S 

0SU10670 

bmk=?.»fmi/x-fm 

OSUI 0630 

IF  IN  11 80, 170,180 

OSUIOAAO 

I 70 

iu=eNK 

OSUI  06  50 

180 

ALPHA:ALPHA48MK 

OSUI 0660 

B J=RJ/ALPHA 

OSUI 0670 

IF (ABS(BJ-HPREVI-APS(0*BJI )200,?00,190 

OSUI06B0 

I»0 

BPRFV=BJ 

OSUI 0690 

IIK-1 

OSUI0500 

200 

RF  TORN 

0SU105I0 

FNP 

OSU10520 

r 

nsul 0530 

c 

. .n^u  10*5  40 

f 

OSUI 0550 

r 

SUBRODTlNf  BFSV 

OSUI 0560 

c. 

OSU10570 

c. 

PURPOSf 

C1SU105B0 

(. 

COMPUTE  THE  V BESSEL  FUNCTION  FOB  A GIVEN  ARGUMENT  AND  ORDFROSUl 0590 

r. 

nsui06P0 

r 

USAGE 

nsuioMO 

f. 

CALL  BFSYI X,N,BY, IFR 1 

OSUI0670 

f 

OSU10630 

r. 

DFSCRIPTION  OF  PARAMETERS 

DSU10660 

c 

X -THF  ARGUMENT  OF  THE  Y BFSSFL  FUNCTION  OFSIRED 

OSUI  0650 

(. 

N -THF  OROFR  OF  THE  Y BESSEL  FUNCTION  DESIRED 

OSUI 0660 

r 

BY  -The  be  suit  ant  Y BESSEL  FUNCTION 

OSUI  0670 

( 

IFR-RE SULTaNT  ERROR  CODE  WHERE 

fiSIM  0^B0 

c 

IFR  = 0 NO  ERROR 

0SIM0640 

c 

lERtl  N IS  NEGATIVE 

OSUI 0700 

f 

IFR=?  X IS  NEGATIVE  OB  TFRO 

OSUI 07 10 

r. 

IER=3  BY  MAS  FXCEEDFD  MAGNITUDE  OF  I0**70 

OSUIO770 

r. 

OSU10730 

c 

RFMARKS 

OSUI 0760 

c 

VERY  SMALL  VALUES  OF  X HAY  CAUSE  THE  RANGE  OF  THF  LIBRARY 

0SU10750 

r 

function  AlOG  TO  BE  EXCEEDED 

OSUI 0760 

l“. 

X MUST  BE  GREATER  THAN  IFRO 

DSUI0770 

f 

N MUST  BE  GREAlFR  THAN  OR  EQUAL  TO  ZERO 

OSUI07B0 

r 

OSUI 0790 

f 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  RFOUIRFD 

OSUlOfiOO 

c 

NONE 

OSUI OB  10 

c. 

nsuiofl?o 

r. 

MF  TMIIO 

nsiiir>pTO 

94 


r%r>r»  r^r-r*  rvno  or*r>0'^r»r»r>r*. 


I'f 


HFCUKHfNCf  RtL*1I0N  *ND  POLYNOMIAL  APPROXIMATION  TFCHNIOIII 
Ai  msCRlflfO  OY  A . J.M.HI TCHCOCK , •POL YNOMI AL  APPROXIMATIONS 
10  BfSSFI  Fl/NC  riDNS  OF  OHOFR  ZFRO  ANO  ONF  AND  TO  RFLATFO 
function;.*,  M.T.A.C.,  V.n  .lOSY.PP.flB-flfi,  ANO  G.N.  WATSON, 
•A  THIATISF  ON  THI  THFORY  OF  RFS<FL  FUNCTIONS',  CAMBKIOr.T 
UNIVFRSlTY  PRFSS,  1RS8,  P.  6? 


SUhROUTINF  BFSYI X,N,BY,IER) 

CHFCX  FOR  FRRORS  IN  N ANO  X 

IF (NllflO, 10, 10 
10  ICR=0 

IF  (XI 190 , 190, ?0 

BRANCH  IF  X LFSS  THAN  OR  EQUAL  A 

?0  IF ( X-A .0 l*O,A0,30 

CUMPUTf  YO  AND  VI  FOR  X GRFATFR  THAN  A 

30  T1=A.0/X 
T2=T1*T1 

P0*( ( ( (-.000O0370A3*T?*.00O0173565)*T2-.0O0OABT613)*T2 
1 ♦.000173a3)*T2-.OOI753062I*T?« .3989A23 

Q0=( I ( ( .0000032312*T2-.00001A2078l»T2A.00003A2AfeB)*T2 
1 -.0000869791 )«T2».OOOAS6A32A IAT2-.01 2A669A 

PI -I ( ( ( .00000A2A 1A*T2-.00002009201*T2*.OOOOS607S9)»T2 
1 0002 23203  I *T ?♦. 00292  I 8 26  I* T 2 *.3989^23 

01*-(  ( ( (-. 0000036 69A*T2*.00001622I»T2-.000039B70BI*T2 
1 ♦.000IO6A7A1 1*T2-.0006390A00I*T2«.037A008A 

A=2.0/S0RTIX) 
h=A»T  I 

C=X-.78S39B2 

YO=A*P0»SIN(CI*B*00*C0S(C) 

Y l=-A*Pl»COS(C I+BPQIASINIC ) 

GO  TO  90 

COMPUTE  YO  ANO  Y1  FOR  X LESS  THAN  OR  EQUAL  TO  A 

AO  XXeX/2. 

X2'XX»XX 

T=ALnC(XX|.. 6772167 
SUM=0. 

Tf  RM=T 
YO-T 

on  70  L*1 ,16 
IF lL-1 160,60,60 
60  SUM»SUM« 1 ./FLOAT (L-1 ) 

60  FL'L 

TS=T-SUM 

TfRM=(TERM*(-X2)/FL**2)*II.-l./(FL*TS)l 
70  VO*YO«TFRM 

TERM  = XX*(T-.6) 

SUM*0. 

YI  » IFRM 
00  BO  L»?,16 

suM*suM« 1 ./float (L-1 I 
FL  ^L 

FL 1»FL-1 . 

T SsT-SUM 

TFRM,^  ( IFRMA  (-X2I/(FL  1*FL  1 )♦(  (TS-.5/FL  )/(TS».6/FLl  I I 
80  Yl=YI«Tf RM 
PI  2».636619fl 
Y0»Pi;*Y0 
YI*-PI2/X«PI2*Y1 


OSUIOBAO 

OSU10B60 

nsui  01,60 
ns  U I 08  70 

OSUl 0880 

OS  III  or,  90 

OSU10900 
OS  U I 09  10 
nsui 0920 
OSUI  09  30 

nsui 09AO 

OSUl 0960 
OSUl 0960 
OSUl  0970 
OSUI0980 
0SU10990 

nsui 1000 
nsui 1010 
nsui 1020 
nsui  1030 
nsui lOAO 
nsui 1060 
OSUl 1060 
nsui 1070 
OSUl  1080 

nsui 1090 
nsui 1100 
nsui  1 1 10 
nsui 1120 
nsui  1130 
OSUl 1 lAO 

nsui 1160 

OSUl 1160 

nsui  1 1 70 
nsui 1180 

OSUl  1190 
OSUl 1200 

nsui 1210 

OSUl 1220 
nsui 1230 
OSUl I 2A0 

nsui  1260 
nsui 1260 
nsui  1270 
OSUl  12H0 
nsui 1290 
nsui 1300 
OSUl 1310 
nsui 1370 
nsui  1330 
OSUl 13A0 
OSUl 1360 
OSUl 1360 
nsui 1370 
nsui  1380 
nsui  1390 

OSUl lAOO 

nsui lA  10 
nsui  I A20 
nsui  1A30 
nsui  lAAO 
nsui  IA60 
nsui IA60 

OSUl  1A70 

nsui 1A80 

liSUl  1A90 

nsui  16(0 

OSUl 1610 
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PC 

100 

no 

no 

1 30 

no 

m 

ns 

ISO 

1 60 
1 70 
1 HO 

190 


CHFCK 

IF  ONLY 

YO 

OR  Y1 

IS 

DFSIRFD 

IF (N-1 

1 100, 100,1 30 

RE  TURN 

EITHER 

YO 

OR  Y1 

AS 

RFOUIREO 

IF (NI 1 
PY=Y1 
GO  TO 
PY=Y0 
GO  TU 

10,170, 

170 

170 

1 10 

orRFOP»3  RECURRFNCF  (lOFRAFIONS  TO  FIND  VN(X) 

YA=YC 
Y0  = YI 
K = 1 

T=FL0AT( 7*K)/X 
YC=I*YB-Y* 

IF  lAHsiYC  )-i.0F7oms,ns,ni 
IF«=3 
RFTIIRN 
K = K»1 

IF  (K-NUSOt  160,1  SO 

YA  = YH 

YP=YC 

r,n  TO  no 

PY  = YC 

RETURN 

IFR=I 

Rt  TURN 

IER  = ? 

RETURN 

FNO 


SUPROUTINF  OSF 
PURPOSE 

TO  compute  tuf  vector  of  integral  values  for  a given 

FQUIUirTANT  TARLE  OF  FUNCTION  VALUES. 

USAGE 

CALL  OSF  IU,Y,Z,NOIM) 

OFSCRIPTION  OF  PARAMFTFRS 
M 
Y 
Z 


NOIM 


THE  INCREMENT  OF  ARGUMENT  VALUES. 

thf  input  vector  of  function  values. 

the  resulting  vector  of  INTFE.RAL  values.  7 MAY  RE 

lOFNTICAL  with  Y. 

The  dimension  of  vectors  y and  7. 


REMARKS 

NO  ACTION  IN  CASE  NOIM  LESS  THAN  3. 

SUBROUTINES  AW  FUNCTION  SUBPROGRAMS  PFOUIRFO 
NONE 

method 

BECINNINC  WITH  7M)*0,  EVALUATION  OF  VFCIOR  7 IS  DONE  BY 
MEANS  OF  SIMPSONS  RULI  TOGETHER  WITH  NEWTONS  3/H  RULE  OR  A 

combination  of  thfsf  two  rulfs.  truncation  error  is  of 

tJRDFR  H*»S  lI.E.  FIXIRTM  OROFP  MFTHOO).  ONLY  IN  CASF  N0IM  = 3 
TRUNCATION  ERROR  OE  7(71  IS  OF  ORDER  H**A. 

FOR  RFFFRFNCF,  SFF 

(II  F .B. HI  IDE  BRAND,  INTRODUCTION  TO  NUMFRICAl  ANALYSIS, 


OSUI |S?0 
OSUI 1530 
OSUI ISAO 
OSUI ISSO 
OSUI IS60 
OSUI IS70 
OSUI ISBO 
OSUI IS90 
OSUI 1600 
OSUI  1610 
OSUI 1670 
OSUI  1630 
OSUI 16*0 
OSUI  I6S0 
OSUI 1660 
OSUI 1670 
OSUI I6B0 
OSUI  1690 
OSUI 1700 
OSUI 1710 
OSUI  1 770 
OSIPl  1730 
OSUI 17*0 
OSUI 1750 
OSUI 1760 
OSUI 1770 
OSUI  1 7P0 
OSUI 1790 
OSUI  IHCO 
OSUI IPIO 
OSUI  1"70 
OSU1IB30 
OSUI  16*0 
OSUI  IB50 
OSUI TH60 
OSUI  1670 
, OSUI 1060 
OSUI IH90 
OSUI 1900 
OSUI  1910 
OSUI 1970 
OSUI n 30 
OSUI 19*0 
OSUI 1950 
OSUI I960 
OSUI  1970 
OSUI loPO 
OSUI 1990 
OSUI  7000 
0SU17010 
OSUI  7070 
OSUI  7070 
OSUI 70*0 
OS U17050 
0SU17060 
OS  II 1 70  70 

nsui706o 
OSUI  70'>0 
OSUI 7100 
OSUI  71  10 
0SU17I70 
OSU17I30 
OSU17i*0 
OSUI  71  50 
OSUI  71  60 
OSUI  71  70 
OSUI 71  BO 
OSUI  7 no 
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oo  rvr»  n r<o  rtn'^n~>nn 


r 

r 


1 


3 


7 


h 


Q 


'<(.rnAW-(MLL,  NfW  vnKK/TO«nNIO/LnNr)nN,  l<y?>6,  PP. 71-76. 
(?)  «./llPMUfHL,  PRAKTISCM)  MATHFMATIK  FllfP  TNCrNIFIIBf  DND 

PHVSIKFR,  ?PB|N(,FR,  BfRLIN/f,OKTTI>4r,fN/HF  inFlPFRG,  1P63, 
I’P.il'.-??!. 


SUBRnUTINF  OSF )H,Y,7,NniHI 


muFNSlON  Y)  1 ) ,Z  (1) 

HT=.333333»*H 
IF (N01M-‘ 17, fl, 1 

Nf)IM  IS  r.RFATFR  THAN  S.  PRFPARATIONS  OF  INTFGRATION  LOOP 
Sl)Ml-YI7)*V(?) 

Sl)Ml  = SllMl  ^SU^^1 
suMi-Ht* ) Y(  n«siim «Y( 3) ) 

AUXl=Yli.)«Y(i.l 
AlIXIiAUXi  «AIIX1 

AIIX1  = SUMI»HT*(Y(3)«AUXI4YI5)I 

A(IX?*HT»  ( Y(  1)4  3.  S7->*(  V(  2 |4Y(5)  I 4 2.6?  5*  ( Y ( 3 ) 4 Y ( <• ) ) 4 V ( 6 1 1 
SUM?=Y ( S ) 4Y( S) 

SUM?  = SllM2  fSIIM? 

SUM2-AUX?-Hl»(  Y(  4Sl)M?4Y(6)  1 

/ ( n = o. 

AUX^YI 31 4Y( 3 ) 

AOX*  AIJX4AI»X 

7 (2  I = MJH?-H7 ♦) Y(?)4AUX4Y(A)) 

Z(31  = SllM) 

?(Al  = Sll>32 
1 F INDlM-6 )S ,6,2 

TNirCPAlION  LOOP 
no  <.  I - 7,  NO  I H,  2 
SOMl=AOXl 
SOM?=AOX2 

AUX1=Y ( I-l 1 4Y(  I-)  1 
AOXI^AOXl 4A0XI 

AIIX1  = S0M1  4MT*(  Y(  l-?)4A0Xl4Vn  1 1 
2(  I-?I  = SUM| 

If  ( l-NOI N )3,6,6 
AOX2=-Y  ( 1 )4Y(  1 1 
AOX?*AUX24AUX? 

AIIX?=SOM2  4HT»(Y(I-I)4AUX24Y)l4l)) 

? I I-l ) =SOM7 
7 )nniM-i  ) = Aoxi 
7 (NOIM 1»AUX2 
Rl  TORN 

/(NOlM-l  MSOM? 

7(ri(llM)  = AUXl 
Rf  TURN 

FNC)  OF  INTFGRATION  LOOP 
I F (NoiM-3  n 2,  n . p 
):niM  IS  rcoAi  lo  a or  s 

SUM2=  I.I?''4Hl*(Yri)4Y(2)4Y(7)4y(2)4Y(3)4YI3)4Yn)4Yt«l) 

SOIU  = Y)  2 ) 4Y(  7 I 

S0M1'SIIM14<  UNI 

S0M!*HTA ( Y(I )4S0M1 4Y( 3) 1 

7 n »40. 

AUXl  =Y  ( 3 )4Y  ( 1) 

AlJXl  4 AUX  1 4 AOX  1 

/( 2)4S0N7-U1 *( Y( 2) 4AOX1 4Y( A) 1 
II  (NIUM-'  nc.p.N 
AUXI  X Y I A ) 4Y  ( <.  ) 


OSOl 2200 

nsol22)0 

OSOl 2220 
nsoi 2230 
nsui?2A0 
nsul 2760 
nsoi  2260 
nsui  2270 
nSM122P0 
nsui 2290 
nsoi 2300 
nsul 2310 
nstii  ?3?o 

0S0123  3P 

nsul 23A0 
nsul 23SO 
nsoi  2360 
nsul 2370 
nsoi  ?3fio 
nsui  2390 
nsoi  ?AOfi 
nsui 2A 10 
nsui  2A  20 
nsoi  2A30 
nsui 2a*o 

OSU12ASO 

nsui  ?A  60 

nsui2A7o 

nsui?A«o 
nsui 2A90 
nsui 2S00 

0SU12S10 
nsui 2S20 
OSUl 2S30 
0SU12SA0 

nsul  2S60 
nsui 2S60 
nsui  2S70 
nsui 2FH0 
nsui 2S90 
nsui 2600 
nsui 26)0 
nsui  26  20 
os  111  26  30 
rsui 26A0 
nsul  2660 
nsul 2660 
nsui  26  70 
OSU126BO 

nsui 2690 
nsui  2700 
OSUl  2710 
nsui  27  20 
nsul 2730 
rsul  27A0 
nsui 27S0 
nsui 2760 
nsui  2 7 70 
nsui 27H0 

IlSlIl  2790 

nsui 2300 
nsui7Hio 
nsoi 2R20 

OSUI2P30 

nsui?PAO 
nsui  2RS0 
nsui ?«60 
nsui  26  70 
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MISSION 

of 

Rcme  Air  Development  Center 


RADC  plaits  and  oondacts  researcb,  exploratory  and  advanced 
development  programs  in  command,  control,  and  oommmicatlons 
(C^)  activities,  and  in  the  areas  of  information  sciences 
and  Intelligence.  The  principal  tectnical  mission  areas 
are  caamunicatioas , electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  scimces,  micramave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


